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engineer  and  technical  monitor  for  the  Naval  Air  Development  Center 
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SUMMARY 


A digital  computer  program  has  been  developed  for  simulating 
the  launch  and  recovery  of  VSTOL  aircraft  operating  from  small 
ships.  The  program,  known  as  VOLAR  (Vought  Launch  and  Recovery 
Dynamics  Program),  employs  the  computational  technique  of  non- 
linear covariance  propagation.  This  permits  the  time  histories 
of  the  means  and  variances  of  all  system  state  variables  to  be 
computed  from  a single  run,  as  opposed  to  the  Monte  Carlo  technique, 
which  requires  multiple  runs,  plus  subsequent  averaging.  Typically, 
VOLAR  requires  approximately  7 percent  of  the  computer  time  required 
for  comparable  Monte  Carlo  simulations. 


The  program  includes  a general  airframe  mathematical  model 
suitable  for  helicopters  or  fixed-wing  aircraft.  The  airwake  of 
a small  ship  is  also  modeled  and  ship  motion  models  are  included 
Two  alternative  models  of  the  human  pilot  are  supplied,  one  is 
based  on  verbal  adjustment  rules,  the  other  is  based  on  optimal 
control  theory,  utilizing  performance  index  parameters  deduced 
from  manned  simulator  experiments.  The  program  is  demonstrated 
for  AV-8A  recovery  on  a small  ship.  The  trends  predicted  by 
VOLAR  are  shown  to  agree  with  flight  test  data. 


Volume  I of  this  report  describes  the  program.  Volume  II  is 
the  user's  manual . 
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A.  1 Philosophy 

VOLAR  is  designed  to  be  a useful  analytical  tool  for  any  task  involving 
covariance  propagation.  It  definitely  is  not  limited  or  restricted  to  the 
examples  presented  in  this  user's  guide.  VOLAR  is  not  a simple  routine  to 
learn,  but  once  mastered  it  provides  the  familiar  user  with  a powerful 
analytical  tool.  The  intended  user  should  have  knowledge  of  differential 
equations,  matrix  algebra,  aerodynamics,  modeling,  control  system  analysis 
and  synthesis,  and  Fortran  programming. 

VOLAR  is  a relatively  flexible  program.  The  basic  structure  allows  one 
to  tailor  the  routine  to  his/her  particular  need  or  requirement.  Then  once 
a problem  is  run,  the  user  may  select  from  several  types  of  output;  including 
line  printer,  printer  plots,  CalComp  plots  and  punch  cards. 

The  program  is  written  in  Fortran  IV  and  COMPASS  for  CDC  6000  or  CYBER 
series  computers.  Each  COMPASS  subroutine  has  a Fortran  equivalent. 

The  subroutines  include  many  comment  statements;  and  for  some  users,  will 
serve  as  their  own  best  user's  guide.  Great  care  was  taken  to  program  in  a 
simple  concise  manner.  When  complex  programming  appears,  the  authors  felt 
it  was  necessary  due  to  speed  and/or  efficiency  considerations. 

A. 2 Some  Comments  on  the  Use  of  VOLAR 

The  user's  guide  should  be  regarded  as  just  that  - a guide.  VOLAR 
cannot  be  completely  learned  from  a book.  The  prospective  user  must  get 
hands-on  experience.  It  is  highly  recommended  that  VOLAR  be  activated  and 
exercised  thoroughly.  The  two  sample  problems  presented  in  this  document 
provide  a good  starting  point. 

The  first  step  in  the  learning  process  is  to  review  the  general  purpose 
and  function  of  each  subroutine  (Section  A. 3).  Next,  the  user  should  read 
Section  A. 4 which  introduces  the  I and  R arrays  and  then  become  familiar 
with  the  main  program,  VOLAR,  and  its  options  (Section  A. 5).  After  this, 
select  one  of  the  sample  cases  and  compare  the  problem-specific  programming 
(refer  to  Appendices  B and  C)  in  the  subroutines  (particularly  DEQU,  PROP, 
and  SETUP)  with  the  general  descriptions  in  Section  A. 3.  The  explanation 
of  the  input  decks  in  Section  A. 8 or  Section  A. 9 is  next.  Following  this. 
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activate  VOLAR  on  your  system.  After  running  the  sample  cases,  try  changing  a 
gain  or  time  constant;  CalComp  plot  a new  variable:  or  change  output  forms. 

The  important  point  is;  roll  up  your  sleeves  and  dig  in. 

A. 3 Description  of  Subroutines 

This  section  presents  a brief  description  of  the  subroutines.  There 
are  basically  two  types  of  subroutines.  The  first  type  falls  into  the  broad 
category  of  Support".  These  subroutines  perform  canned  operations  (e.g.  matrix 
multiplication,  integration,  etc.)  and  do  not  change  from  problem  to  problem. 
The  second  type  of  subroutine  is  very  much  problem  dependent.  Through  these 
subroutines  the  user  tailors  VOLAR  to  model  and  solve  his  particular  problem. 
Subroutines  requiring  user  interface  are  so  noted  in  the  descriptions  which 
follow. 

ADAMS  - This  subroutine  performs  second  order  Adams 
integration.  It  is  the  simplest  of  the  three 
integration  techniques  available  to  the  user. 

All  the  integration  subroutines  are  constructed 
to  integrate  a set  of  first  order  differential 
equations  describing  the  means  and  three 
covariance  matrices  (three  are  required  when 
an  optimal  pilot  is  used).  The  calling 
routine  is  VOLAR. 

AIRWAKE  - This  subroutine  calculates  airwake  disturbance, 
means  and  standard  deviations,  due  to  the  ship's 
presence.  The  programmed  eauations  are  from 
Reference  12.  The  calling  routine  is  DEQU. 

AIRWT8L  - The  airwake  data  tables  of  Reference  12  are 

contained  here  in  DATA  statements.  This  subroutine 
is  called  once  by  SETUP  to  construct  a case-particular 
particular  (VWqD,  *wod)  set  a^rwake  tables  to 
be  used  by  AIRWAKE. 
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AXIS 

CALPLT 

DASPLT 


This  subroutine  produces  a coordinate  axis  for 
CalComp  plots.  The  calling  routine  is  CALPLT. 

This  is  a general  purpose  CalComp  plot  routine 
The  calling  program  is  VOLAR. 

This  subroutine  draws  dashed  lines  on  CalComp 
plots.  The  calling  routine  is  CALPLT. 


DEQU 


DESCRIB  - 

I 


ELLIPSE  - 


This  subroutine  contains  the  differential 
equations  describing  the  propagation  of  the 
means,  the  aero  data  table  look  ups  and  open 
loop  control  inputs.  The  reference  flight  values 
about  which  the  trajectory  is  linearized  are 
defined  and  updated  here.  DEQU  is  one  of  the 
user  defined  subroutines.  It  must  be  programmed 
by  the  user  to  meet  the  requirements  of  his 
particular  application.  DEQU  is  called  by  SETUP 
and  the  integration  subroutines  ADAMS,  RUNGE, 
and  KUTTA. 

This  routine  calculates  single-valued  nonlinear 
describing  functions  of  one  variable.  TABRD 
must  be  called  before  the  first  call  to  DESCRIB. 
The  calling  routine  for  DESCRIB  is  DEOU. 

This  routine  draws  a 50  percent  probability 
ellipse  on  CalComp  trajectories.  These  ellipses 
can  be  scaled  differently  from  the  trajectory 
to  increase  the  clarity  of  presentation.  The 
calling  routine  is  CALPLT. 

This  routine  evaluates  the  error  function 
associated  with  the  normal  distribution  curve. 

The  calling  routine  is  DESCRIB. 

One  of  three  integration  subroutines.  This 
subroutine  performs  third  order  Runge-Kutta 
integration.  The  calling  program  is  VOLAR. 
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MATINV 


MATRIX 


MISCAL 


PRMPLT 


ROTATE 


This  subroutine  calculates  the  inverse  of  a 
matrix.  The  calling  routine  is  SOLVR. 

This  routine  handles  general  matrix  operations. 

The  operations  include  adding  two  matrices, 
subtracting  two  matrices,  multiplying  two 
matrices,  multiplying  one  matrix  by  the  transpose 
of  another  matrix,  matrix  transposition,  and  laying 
a matrix  into  another  matrix.  MATRIX  is  programmed 
in  COMPASS.  The  FORTRAN  equivalent  operations 
appear  as  comments  in  MATRIX.  The  COMPASS  version 
executes  about  20  percent  faster  than  the  FORTRAN 
version.  The  calling  routines  are  PROP  and  SOLVR. 

This  subroutine  handles  miscellaneous  calculations 
that  have  no  logical  place  to  go.  MISCAL  is  called 
by  VOLAR  before  the  output  file  is  written.  This 
subroutine  is  user  specified. 

This  is  a general  purpose  line-printer  plot  routine. 
A 7 by  7 inch  plot  is  produced.  The  calling  routine 
is  VOLAR. 

This  subroutine  contains  the  matrix  differential 
equation  for  covariance  propagation.  The  user 
must  define  the  elements  of  the  F and  G matrices 
in  PROP.  These  matrices  define  the  linearized 
system  equations  and  are  normally  composed  of 
several  smaller  submatrices.  PROP  is  called  by 
SETUP  and  the  integration  subroutines  ADAMS, 

RUNGE,  and  KUTTA. 

This  subroutine  rotates  a coordinate  system 
through  any  angle.  The  calling  routines  are  CALPLT 
and  ELLIPSE. 


RUNGE 


SETUP 


SOLVR 


TABRD 


- This  subroutine  performs  fourth  order  Runge- 
Kutta  integration.  It  is  the  more  accurate 

but  most  time-consuming  of  the  three  integration 
methods  available  in  VOLAR.  The  calling  program 
is  VOLAR. 

- This  subroutine  initializes  the  integrators. 

Any  calculations  to  be  initialized  and  one-time- 
only  calculations  also  appear  in  SETUP.  For 
example,  ship  scaling  factors  and  wind  over  deck 
are  calculated  here  because  they  are  fixed  for 

a given  run.  Because  of  its  problem  dependence, 
SETUP  is  user  defined.  SETUP  is  called  by 
VOLAR. 

- This  subroutine  solves  for  the  steady-state 
(t-+«)  solution  of  the  matrix  Riccati  eouation. 

It  may  be  used  to  initialize  the  covariance 
matrix,  P,  before  running  a time  history.  It 
will  not  work  when  more  than  one  covariance  matrix 
(e.g.  the  optimal  pilot  model  case)  is  involved. 
The  calling  routine  is  PROP. 

- This  subroutine  decomposes  nonlinear  functions 
into  even  and  odd  describing  function  elements. 

The  inputs  to  TABRD  are  a table  of  independent 
variables,  a table  of  dependent  variables,  and 
the  number  of  variables  in  the  tables.  TABRD 
must  be  called  once  for  each  nonlinear  element 
before  DESCRIB  can  calculate  the  appropriate 
describing  function.  The  calling  routine  is 
SETUP. 
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TLU  - This  routine  is  used  to  look  up  a series  of 

dependent  variables.  Linear  interpolation  is  used 
to  generate  the  dependent  variables.  With  one 
call,  several  dependent  variables  can  be  made 
available  to  COMMON  provided  that  all  the 
dependent  variable  tables  are  associated  with 
the  same  independent  variable  table.  The  calling 
subroutines  are  DEOU  and  PROP. 

VOLAR  - This  routine  is  the  main  program  and  handles 

almost  all  the  program  input  and  output  through 
its  various  options. 

A. 4 Description  of  the  I and  R Arrays 

The  I and  R arrays  contain  all  the  integer  and  real  variables  of  the 
program  respectively.  The  I array  is  equivalenced  to  the  elements  of  the 
labeled  coimion,  "INTEGER".  The  R array  is  equivalenced  to  elements  of 
blank  common.  The  arrays  not  only  handle  the  communication  between  the 
subroutines  but  also  all  of  the  programs  input  and  output.  The  I and  R 
arrays  are  user  defined. 

A. 4.1  Program  Modification 

If  the  user  wishes  to  modify  the  program,  any  new  variables  would  be 
added  to  the  common  blocks  (I  or  R arrays).  This  gives  the  user  the  capability 
of  inputting  or  outputting  the  variables.  The  new  variables  may  be  added  any 
where  within  the  commons. 

A. 4. 2 Input  and  Output  through  the  I and  R Arrays 

All  of  the  programs  input  and  output  were  designed  around  the  I and  R 
arrays.  Since  the  I and  R arrays  (common  blocks)  have  many  variables  and 
arrays  within  them,  bookkeeping  the  variables  is  difficult.  To  simply  the 
bookkeeping,  the  I and  R arrays  are  divided  into  segments  and  elements.  Each 
segment  Is  defined  by  a specific  number  of  elements.  The  segments  are 
classified  by  groups  of  elements  which  are  defined  by  a collection  of  simple 
variables  or  subscripted  variables.  The  lengths  of  the  segments  are  specified 
in  the  problem  initialization  phase  (A. 5.1).  The  initialization  phase  actually 


A- 7 


aligns  the  commons  with  the  I and  R arrays.  Care  must  be  taken  in  setting 
the  lengths  of  the  segments  after  program  modification  because  the  input  and 
output  will  be  malaligned.  Figure  A.l  is  an  example  of  a simple  R array.  In 
this  example,  DTOUT  may  be  referred  to  as  segment  1,  element  4.  The  fifth 
element  of  TABV2  is  segment  3,  element  5.  Once  the  user  has  established  the 
order  of  variables  in  the  R array,  then  each  variable  may  be  referred  to  by 
its  segment  number  and  element  number  within  that  segment.  This  facilitates 
input/output  and  simplifies  modifications  to  the  program.  Understanding  the 
concept  of  the  I and  R arrays  is  extremely  important  to  the  use  of  VOLAR. 

C0MM0N  TIME,  T0,  DTIME,  DT0UT , DT0UT2,  TMAX 

C0MM0N  TABV1(6),  TABV2(10) 


Variable 

TIME 

T0 

DTIME 

DT0UT 

DT0UT2 
TMAX 
TAB  VI 
TABV2 


Segment  No. 
1 
1 
1 
1 
1 
1 
2 
3 


Element  No. 
1 
2 

3 

4 

5 

6 

1-6 

1-10 
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Figure  A.l 


li 


Appendix  E presents  the  variables  in  the  R and  I arrays  for  the  analysis 
reported  in  this  document. 
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A. 5 The  Main  Program  (Description  of  Options) 

The  program  is  divided  into  two  input  phases.  The  first  phase  is  program 
initialization  which  defines  the  lengths  for  the  R and  I arrays.  The  second 
phase  is  user  defined  and  directs  the  program  through  the  various  program 
options . 

A. 5.1  Program  Initialization  (Phase  I) 


Definition  of  the  R and  I array  lengths  are  input  through  the  following 
namelist.  The  lengths  of  the  segments,  defined  in  the  program  input,  must  be 
compatible  with  the  number  of  variables  in  blank  and  labeled  COMMON.  The 
segment  lengths  allow  the  main  program  to  bookkeep  the  absolute  location  of 
a variable  in  COMMON. 


A. 5. 2 Program  Options  (Phase  II) 

The  program  is  divided  into  28  options  to  handle  problem  execution,  input, 
and  output.  The  definition  of  the  options  are  as  follows: 

Option  Task 

1 Write  I and  R array  I/O  instructions  on  a file. 

2 Store  data  into  the  I array  and  write  I array  data  on  a file. 

3 Print  data  from  the  I array. 

4 Store  data  into  the  R array  and  write  R array  data  on  a file. 

5 Print  data  from  the  R array. 

6 Setup  the  differential  equations. 

7 Designate  the  time  history  variables  from  the  R array  to  be 

written  to  a file  by  option  21. 

8-9  Inactive. 

10  Read  data  from  input  and  generate  a file  compatible  for 

options  14,  15,  17,  or  18. 

11  Store  tabulating  or  punching  instructions  on  a file.  The  file 
will  be  used  in  option  14  or  15. 

12  Store  printer  plotter  instructions  on  a file.  The  file  will  be 
used  by  option  17. 

13  Store  Cal  Comp  plotter  instructions  on  a file.  The  file  will  be 
used  by  option  18. 

14  Print  the  time  history  output  in  a tabular  form. 

15  Punch  cards  from  the  time  history  output. 

16  Inactive. 

17  Plot  the  time  history  output  on  the  printer. 

18  Plot  the  time  history  output  on  a CalComp  plotter. 

19  Print  the  contents  of  a file  in  blocked  form. 

20  Inactive. 

21  Generate  a time  history. 

22-27  Inactive. 

28  Store  the  R array  on  Tape  28  or  retrieve  the  R array  from 
Tape  28. 


Input  the  following  card  if  I < 0. 


4. 

1-40 

41-80 

FMTI 

FMT0 

where: 


NOPT 

- Option  number  (input  a 1 for  this  option). 

(15) 

LBL0UT 

- Description  of  option. 

(7A10) 

NLIS7 

- File  where  the  instructions  will  be  written. 

(15) 

I 

- Beginning  segment  number  of  the  I or  R array. 

(15) 

J 

- Beginning  element  number  of  the  I or  R array. 

(15) 

K 

- Ending  segment  number  of  the  I or  R array. 

(15) 

L 

- Ending  element  number  of  the  I or  R array. 

(15) 

LABEL 

- Variable  names  for  the  elements  defined  between 

1 ,J  and  K,L. 

(5A10) 

SCALE 

- Scale  factor  to  be  applied  to  the  elements. 
Default  is  one. 

(E10.3) 

FMTI 

- Variable  input  format. 

(4A10) 

FMT0 

- Variable  output  format 

(4A10) 

A- 11 


NOTE: 


1.  Multiple  sets  of  cards  3 and  4 may  be  input. 

2.  FMTI  and  FMT0  are  input  only  if  I < 0. 

3.  Option  1 is  terminated  when  a blank  card  is  input  for  I,J,K,L. 

4.  See  Paragraph  A. 5. 3. 2 for  a sample  input. 

A. 5.2.2  Option  2 - Store  data  into  the  I array 


Input  the  following  cards  if  NLIST  = 5. 


3. 

1-5 

6-10 

11-15 

16-20  21-70 

71-80 

I 

J 

K 

LABEL 

SCALE 

Input  the  following  card  if  I < 0. 


where: 

N0PT  - Option  number  (input  a 2 for  this  option).  (15) 

LBL0UT  - Option  description.  (7A10) 

NLIST  - File  containing  input  instructions.  (15) 


0 

D 

li 

[] 


Li 

” 

•• 

I 

I 

n 

n 

o 

n 

a 

o 
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NDATA  - File  containing  I array  data. 

NPF  - File  to  be  created  containing  I array  data. 

If  NPF  = 0,  no  file  is  created. 

I - Beginning  segment  number  of  the  I array. 

J - Beginning  element  number  of  the  I array. 

K - Ending  segment  number  of  the  I array. 

L - Ending  element  number  of  the  I array. 

LABEL  - Variable  names  for  the  elements  defined 
between  I,J  and  K,L. 

SCALE  - Scale  factor  to  be  applied  to  the  elements 
in  the  I array  between  I,J  and  K,L. 

FMTI  - Variable  input  format. 

FMT0  - Variable  output  format. 

IN  - 1 array  elements  (Default  is  1615  and  may 

be  changed  to  FMTI). 

NOTE: 

1.  Multiple  sets  of  cards  3,  4 and  5 may  be  input. 

2.  The  format  for  card  5 may  be  changed  by  specifying 
a negative  I on  card  3 and  inputtim  the  format  on 
card  4. 

3.  If  NLIST  and  NDATA  are  not  equal  to  5,  cards  3,  4,  and  5 
are  input  from  files  NLIST  and  NDATA. 

4.  Option  2 is  terminated  when  a blank  card  is  input  for 
1 ,J,K,L. 

5.  See  Paragraph  A. 5. 3.3  for  a sample  input. 


y 

n 
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A. 5. 2. 3 Option  3 - Print  from  the  I array 


1. 

1-5 

6-10 

11-80 

N0PT 

LBL0UT 

2. 

1-5 

NLIST 

Input  the  following  cards  if  NLIST  = 5. 


1-5 

6-10 

11-15 

16-20 

21-70 

71-80 

I 

J 

K 

L 

LABEL 

SCALE 

Input  the  following  card  if  I < 0. 


4. 

1-40 

FMT0 

where: 


N0PT 

- Option  number  (input  a 3 for  this  option). 

(15) 

LBL0UT 

- Option  description. 

(7A10) 

NLIST 

- File  where  print  instructions  reside. 

(15) 

I 

- Beginning  segment  number  of  the  I array. 

(15) 

J 

- Beginning  element  number  of  the  I array. 

(15) 

Y. 

- Ending  segment  rumber  of  the  I array. 

(15) 

L 

- Ending  element  number  of  the  I array. 

(15) 

LABEL 

- Variable  names  for  the  elements  defined 
between  I,J  and  K,L. 

(5A10) 

SCALE 

- Scale  factor  to  be  applied  to  the  I array 
elements  between  I,J  and  K,L  when  printed. 
Default  is  one. 

(E10.3) 

FMT0 

- Variable  output  format. 

(4A1C) 

!J 

D 

D 

D 


l J 


i 

I 

T 

■ 

fl 
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1.  Multiple  sets  of  cards  3 and  4 may  be  input. 


2.  The  print  format  may  be  changed  by  specifying  a 
negative  I and  inputting  the  format  on  card  4. 
The  format  is  maintained  until  it  is  reset  or 
another  option  is  specified. 

3.  The  option  is  terminated  when  a blank  card  is 
input  for 

4.  See  Paragraph  A. 5. 3. 4 for  a samole  input. 


A. 5. 2. 4 Option  A - store  data  into  the  R array,  and  write  the  R 
array  on  to  a file 


1-5  6-10 

N0PT 


LBL0UT 


11-80 


2.  :.-5  6-10  11-15  

NLIST  NDATA  I NPF  ! 

Input  the  following  2 cards  if  NLIST  = 5. 


-5  I 6-10  1 11-15  16-20  2 


BF 


71-80 


D 


Input  the  following  card  if  NDATA  = 5. 

5.  1-10  11-20 
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where: 


N0PT 

- 

Option  number  (input  a 4 for  this  option). 

(15) 

LBL0UT 

- 

Option  description. 

(7A10) 

NLIST 

- 

File  containing  input  instructions. 

(15) 

NDATA 

- 

File  containing  R array  data. 

(15) 

NPF 

- 

File  containing  R array  data  is  to  be  created. 

If  N^F  = 0,  no  file  is  created. 

(15) 

I 

- 

Beginning  segment  number  of  the  R array. 

(15) 

J 

- 

Beginning  element  number  of  the  R array. 

(15) 

K 

- 

Ending  segment  number  of  the  R array. 

(15) 

L 

- 

Ending  element  number  of  the  R array. 

(15) 

LABEL 

- 

Variable  names  for  the  elements  defined 
between  I ,J  and  K,L. 

( 5A 10 ) 

SCALE 

- 

Scale  factor  to  be  applied  to  the  D array 
elements  between  I,J  and  K,L.  Default 
is  one. 

(E10.3) 

FMTI 

- 

Variable  input  format. 

(4A10) 

FMT0 

- 

Variable  output  format. 

(4A10) 

R 

- 

R array  elements  (Default  is  8F10.3  and  may 
be  changed  to  FMTI). 

NOTE: 

1.  The  format  for  card  5 may  be  changed  by  specifying  a negative 
I on  card  3,  and  inputting  the  format  on  card  4.  The  format 
is  maintained  until  it  is  reset  or  another  option  is 
specified. 

2.  If  NLIST  and  NDATA  are  not  eaual  to  5,  cards  3,  4,  and  5 
are  input  from  files  NLIST  and  NDATA. 

3.  Multiple  sets  of  cards  3,  4,  and  5 may  be  Input. 

4.  Option  4 is  terminated  when  a blank  card  is  input  for 
I ,J,K,L. 

5.  See  Paragraph  A. 5. 3. 5 for  a sample  input. 
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Option  5 - Print  from  the  R array. 


■ 


A. 5. 2. 5 


1. 

1-5 

6-10 

11-80 

N0PT 

LBL0UT 

2. 

1-5 

NTAP 

Input  the  following  cards  if  NTAP  = 5. 


1-5 

6-10 

11-15 

16-20 

21-70 

71-80 

I 

J 

K 

L 

LABEL 

SCALE 

Input  the  following  card  if  I < 0. 


4. 

1-40 

FMT0 

where: 

N0PT  - Option  number  (input  a 5 for  this  option).  (15) 

LBL0UT  - Option  description.  (7A10) 

NTAP  - File  where  write  instructions  reside.  (15) 

I - Beginning  segment  number  of  the  P.  array.  (15) 

J - Beginning  element  number  of  the  R array.  (15) 

K - Ending  segment  number  of  the  R array.  (15) 

L - Ending  element  number  of  the  R array.  (15) 

LABEL  - Variable  names  for  the  elements  defined  (5A10) 

between  I,J  and  K,L. 

SCALE  - Scale  factor  to  be  aoplied  to  the  R (E10.3) 

array  elements  when  printed.  Default 
is  one. 

FMT0  - Variable  output  format.  (4A10) 
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NOTE: 

1.  If  NTAP  is  not  eoual  to  5,  cards  3 and  4 are  input  from 
file  NTAP. 

2.  The  default  print  format  of  (7F14.5)  may  be  changed  to 
FMT0  by  specifying  a negative  I and  inputting  the  new 
format.  The  new  format,  remains  in  effect  until  reset  or 
a new  option  is  input. 

3.  Multiple  sets  of  cards  3 and  4 nay  be  input. 

4.  Option  5 is  terminated  by  inputting  a blank  card 
for  I,J,K,L. 

5.  See  Paragraph  A. 5. 3. 6 for  a sample  input. 

A. 5. 2. 6 Option  6 - Set  up  differential  eauations. 

This  option  is  composed  of  five  groups  of  data: 

1.  Derivative  indices. 

2.  State  indices. 

3.  PD0T  and  P indices. 

4.  PED0T  and  PE  indices. 

5.  PXD0T  and  PX  indices. 

Groups  3,  4,  and  5 above  must  be  input  in  the  order  specified.  Each  group 
of  data  is  composed  of  cards  2 and  3 below. 


1. 

1-5 

6-10 

11-80 

N0PT 

LBL0UT 

2. 

1-5 

6-10 

M 

N 

3. 

1-5 

6-10 

11-15 

16-20 

21-70 

I 

J 

K 

L 

LABEL 

il 


3 


1 1 


li 


.. 


D 

D 

D 

2 
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where: 


N0PT 

- Option  number  (inout  a 6 for  this  ontion). 

(15) 

LBL0UT 

- Option  description. 

(7A10) 

M 

- The  segment  number  of  the  I array  where  the 

R array  location  of  the  differential  eouation 
elements  will  be  stored. 

(15) 

N 

- The  element  number  of  the  I array  where  the  R 
array  location  of  the  differential  equation 
elements  will  be  stored. 

(15) 

I 

- Beginning  segment  number  of  the  R array  where 
the  differential  eouation  elements  will  be 
stored. 

(15) 

,1 

- Beginning  element  number  of  the  R array  where 
the  differential  equation  elements  will  be 
stored. 

(15) 

K 

- Ending  segment  number  of  the  R array  where  the 
differential  equation  elements  will  be  stored. 

(15) 

L 

- Ending  segment  number  of  the  R array  where  the 
differential  equation  elements  will  be  stored. 

(15) 

LABEL 

- Variable  names  for  the  elements  defined  between 

I,J  and  K,L. 

(5A10) 

NOTE: 

1.  Multiple  sets  of  card  3 may  be  specified. 

2.  Each  group  of  indices  will  be  terminated  when  a blank 
card  is  specified  for  I,J,K,L. 

3.  PD0T  and  P indices  are  defined  by  I,J  and  K,L  respectively. 

4.  PED0T  and  PE  indices  are  defined  by  I,J  and  K,L 
respectively. 

5.  PXD0T  and  PX  indices  are  defined  by  I,J  and  K,l  respecti vely . 

6.  See  Paragraph  A. 5. 3. 7 for  a sample  input. 
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LJ 


A. 5. 2. 7 Option  7 - Designate  the  time  history  variables  to  be 
output 


1. 

1-5 

6-10 

11-80 

N0PT 

LBL0UT 

2. 

1-5 

6-10 

11-15 

16-20 

21-70 

I 

J 

K 

L 

LABEL 

where: 


ii 

D 

D 

D 


N0PT 

- Option  number  (input  a 7 for  this  option). 

(15) 

LBL0UT 

- Option  description. 

(7A10) 

I 

- Beginning  segment  number  of  the  R array  to 
be  output. 

(15) 

J 

- Beginning  element  number  of  the  R array  to 
be  output. 

(15) 

K 

- Ending  segment  number  of  the  R array  to  be 
output. 

(15) 

L 

- Ending  element  number  of  the  R array  to  be 
output. 

(15) 

LABEL 

- Variable  names  for  the  elements  defined 
between  !,J  and  K,L. 

(5A10) 

NOTE: 

1.  Multiple  sets  of  card  2 may  be  input 

2.  Option  7 is  terminated  when  a blank  card  for 
I,J,K,L  is  input. 

3.  See  Paragraph  A. 5. 3. 8 for  a sample  input. 


D 

D 
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A. 5. 2. 8 Option  10  - Read  data  from  input  and  generate  a file  compatible 
for  options  14,  15,  17,  and  18. 


1. 

1-5 

6-10 

11-80 

N0PT 

LBL0UT 

2. 

1-5 

NLIST 

3. 

1-5 

6-10 

11-15 

K1 

K2 

L 

Input  the  following  card  K1  times. 


4. 

1-10 

11-20 

21-30 

71-80 

R0W(1) 

R0W(2) 

R0W ( 3 ) 

• • • 

R0W ( K2 ) 

where: 

N0PT 

LBL0UT 

Nil  ST 

K1 

K2 

L 


R0W 


Option  number  (input  a 10  for  the  option).  (15) 

Option  description.  (7A10) 

File  where  data  is  to  be  written.  (15) 

Number  of  rows  of  data.  (15) 

Number  of  columns  of  data.  (15) 

Data  Flag  (15) 


L = 0 - Number  of  data  points  eaual  to  the 
number  of  rows. 

L t 0 - First  data  point  indicates  the  number 
of  data  points  in  a floating  point 
number. 

Data  point  array. 


NOTE: 


1.  See  Paragraph  A. 5. 3. 9 for  a sample  inout. 
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0 

Q 
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NOTE: 

1.  K2  is  the  seauential  number  for  the  variable  as  ordered 
in  option  7. 

2.  For  tabulating  instructions,  multiple  sets  of  cards  3 and 
4 may  be  input.  The  option  is  terminated  when  a blank 
card  is  input  for  Kl. 

3.  For  punching  instructions,  a multiple  set  of  card  3 may 
be  input.  The  option  is  terminated  wher  a blank  card  for 
NTAP  and  K2  is  innut. 

4.  The  file  NLIST  will  be  used  in  option  14  or  15. 

5.  See  Paragraph  A.E.3.10  for  a sample  innut. 


where: 


u 

N0PT 

- 

Option  number  (input  a 12  for  this  option) 

• 

05) 

fl 

LBL0UT 

- 

Option  description. 

(7A10) 

1 Li 

NLIST 

- 

File  where  instructions  will  be  written. 

(15) 

11 

1 LJ 

K1 

- 

Number  of  plots. 

05) 

TITLE 

- 

Three  card  plot  title. 

(8A10) 

NCURVS 

- 

Number  of  curves  oer  plot. 

(15) 

XPL07 

- 

Number  of  inches  from  left  paper  margin  to 
left  end  of  the  X axis . 

the 

(F10.3) 

YPL0T 

- 

Number  of  inches  from  top  of  the  paper  to 
bottom  of  the  Y axis. 

the 

(F10.3) 

X 

LU 

- 

Increment  per  inch  or  the  X axis. 

(F10.3) 

I n 

OELY 

- 

Increment  per  inch  on  the  Y axis. 

(F10.3) 

STARTX 

- 

X value  at  axis  intersection. 

(F1C.3) 

STARTY 

- 

Y value  at  axis  intersection. 

(F10.3) 

FMTI 

- 

Label  on  X axis. 

(4A10) 

FKT0 

- 

Label  on  Y axis. 

(4.410) 

Li 

LABEL 

- 

Legend  associated  with  a curve. 

(4410) 

NDATA 

- 

File  where  data  for  the  curve  is  located. 

(15) 

0 

NC0LX 

- 

Column  number  on  the  file  where  tue  Y data 
located. 

is 

(15) 

jo 

NC0LY 

- 

Column  number  on  the  file  where  the  Y data 
located. 

is 

(15) 

n 

NSYM 

- 

Plotting  symbol. 

(Al) 

SCALEX 

- 

Scale  factor  to  be  applied  to  the  X data 
(Default  = 1.0). 

(F10.3) 

SCALEY 

- 

Scale  factor  to  be  applied  to  the  Y data 
(Default  = 1.0). 

(F10.3) 

NOTE: 


1.  See  Paragraph  A. 5. 3. 11  for  a sample  Input. 
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2. 


3. 


4. 


1-5 

K1 


6-10 


11-15 

L 


Input  the  following  cards  K1  times. 


5. 

1 

2-80 

SPL0TD  NCURVS=  ,ANGLE=  ,XPL0T= 

6. 

□ 

2-80 

1 

YP10T=  , XEXIT=  , YEXIT=  , XLEGND= 

7. 

□ 

2-80 

1 

YLEGND=  , X0RG=  , Y0RG=  , LNGTHX= 

8. 

H 

2-80 

LNGTHY=  , STARTX=  , STARTY= , DELX= 


9. 


Input  the  fol 


li 


Input  the  following  cards  NCURVS  times. 


N0PT 

- Option  number  (input  a 13  for  this  option). 

(15) 

LBL0UT 

- Option  description. 

(7A10)  • 

PC0DE 

- Cal  Comp  plot  code. 

(A10) 

11 

NTAP 

- File  produced  by  Cal  Comp. 

(A10) 

jp 

NLIST 

- File  where  plotting  instructions  will  be 
written. 

(15) 

K1 

- Number  of  plots  to  be  produced. 

(15) 

L 

NCURVS 

- Indicator  for  plotting  on  centimeter  paper. 

L = "C.M. " 

- Number  of  curves  per  plot 

(A4) 

0 

D 

O' 
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0 


ANGLE 


Plot  rotation  angle  in  degrees  (0  or  90  degrees). 


J 


0 

0 


0 

0 

n 


XPL0T 

YPL0T 

XEXIT 

YEXIT 

XLEGND 

YLEGND 

X0RG 

Y0RG 

LNGTHX 

LNGTHY 

STARTX 

STAPTY 

DELX 

DELY 

IS0PT 


TITLE 

LARELX 

LABELY 

LEGEND 


Distance  in  inches  for  which  the  origin  of  the 
plot  is  to  be  moved  in  the  X direction. 

Distance  in  inches  for  which  the  origin  of  the 
plot  is  to  be  moved  in  the  Y direction. 

Final  pen  movement  in  the  X direction. 

Final  pen  movement  in  the  Y direction. 

X coordinate  of  the  plot  legend  in  inches. 

Y coordinate  of  the  plot  legend  in  inches. 

X value  at  which  axes  intercept. 

Y value  at  which  axes  intercept. 

Length  of  X axis  in  inches. 

Length  of  Y axis  in  inches. 

Initial  value  on  the  X axis. 

Initial  value  on  the  Y axis. 

Increment  per  inch  on  the  X axis. 

Increment  per  inch  on  the  Y axis. 

Bit  array  for  read  option.  See  Figure  A. 4. 

BIT(C)  = 1 - Read  title  to  be  centered  over  plot. 

BIT(l)  = 1 - Read  X axis  label. 

BIT(2)  = 1 - Read  Y axis  label. 

BIT(3)  = 1 - Read  legend  title. 

BI T (4 ) =1-  Read  coordinates  of  the  title  and 
the  title. 


Four  card  title.  (a5 10) 
X axis  label.  (4A10) 
Y axis  label.  (4A10) 
Legend  for  the  curve.  (4A10) 
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X 


X coordinate  for  the  title. 


(F10.3) 

(F10.3) 


r 


Y 

NDATA 

NC0LX 

NC0LY 

NCX 

NCY 

LINTYP 


DASH 

NSYM 
SCALE X 
SCALEY 
SCALXX 

SCALYY 

K 


Y coordinate  for  the  title. 

File  where  data  to  be  plotted  is  located. 

Column  number  on  the  file  where  the  X data  is 
located. 

Column  number  on  the  file  where  the  Y data  is 
located. 

Column  number  where  the  X elliptical  axis  data 
is  located. 

Column  number  where  the  Y elliptical  axis  data 
is  located. 

Line  type  for  plot. 

LINTYP  = Zero,  line  only  (no  symbol) 

LINTYP  = Negative,  symbol  only  (no  line) 

LINTYP  = Positive,  both  line  and  symbol 

Length  of  plotting  segments  for  dashed  line  in 
inches.  See  Table  A.l. 

CalComp  symbol  code  (99  if  ellipse). 

Scale  factor  for  the  X data  (Default  = 1.0). 

Scale  factor  for  the  Y data  (Default  = 1.0). 

Scale  factor  for  the  X elliptical  axis 

(Default  = SCALEX). 

Scale  factor  for  the  Y elliptical  axis 
(Default  = SCALEY). 

Line  style.  See  Table  A.l. 


LABEL  - Legend  for  line  style. 


NOTE: 

1.  See  Paragraph  A. 5.3. 12  for  a sample  input. 


f | 

i i 


n 


O I 

f 

M i 

u I 

(06) 

(4A10) 


2.  See  Figure  A. 3 for  a plot  of  ellipses  of  probability. 

3.  See  Figure  A. 2 for  a diagram  of  plotting  variables. 
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Table  A.l 


U 


0 

Ln 


LINE 


lstyle 

DASH 

177610 

.01 

177610 

.05 

177430 

.01 

177430 

.05 

177070 

.01 

177070 

.05 

177774 

• 01 

177774 

.05 

177400 

.01 

177400 

.05 

177752 

.01 

177752 

.05 

177714 

.01 

177714 

.05 

125452 

• 01 

125452 

.05 

170360 

.01 

170360 

.05 

177744 

.01 

177744 

.05 

177444 

.01 

177444 

.05 

146314 

.01 

146314 

• 05 

176314 

• 01 

176314 

• 05 

177777 

ANY 
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A. 5. 2. 12  Option  14  - Print  the  time  history  in  a tabular  form 


where: 


N0PT 

- Option  number  (inout  a 14  for  this  option). 

(15) 

LBL0UT 

- Option  description. 

(7A10) 

NLIST 

- File  where  tabulating  instructions  are 
located. 

(IF) 

NOTE: 

1.  File  NLIST  is  generated  ly  option  11. 

2.  See  Paragraph  A. 5. 3. 13  for  a sample  input. 


A. 5. 2. 13  Option  15  - Punch  the  time  history 


where: 


N0PT 

- Option  number  (input  a 15  for  this  option). 

(IB) 

LBL0UT 

- Option  description. 

(7A10) 

NLIST 

- File  where  the  punching  instructions  are 
located 

(15) 
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NOTE: 


1.  File  NLIST  is  generated  by  Option  11. 

2.  See  Paragraph  A. 5. 3. 14  for  a sample  input. 


A. 5. 2. 14  Option  17  - Produce  a printer  plot  of  the  time  history 

1. 


Input  the  following  cards  if  NLIST  = 5. 

3. 


1-5 


| NPL0TS | 

Input  the  following  cards  NPL0TS  times. 


1-5 

6-10 

11-80 

N0PT 

LBL0UT 

1-5 

6-10 

NLIST 

NPRNT 

0 

0 

D 

0 

D 


Input  the  following  card  NCURVS  times. 


1-40  41-45  46-50  51-55  56-57  58  59-60  61-70  71-8C 

LEGEND  NDATA  NC0LX  NC0LY  ISYM  SCALEX  SCALEV 


v/here: 


NCURVS 


NDATA 


- Option  number  (input  a 17  for  this  option). 


LBL0UT  - Option  description. 

NLIST  - File  where  the  plotting  instructions  are 
located. 

NPRNT  - Print  option  to  display  input. 

NPRNT  =0  - no  print 
•NPRNT  = 1 - print 

NPL0TS  - Number  of  plots  to  be  produced. 

TITLE  - Three  card  plot  title. 


- Number  of  curves  per  plot. 

- Number  of  inches  from  the  left  paper  margin 
to  begin  the  X axis.' 

- Number  of  inches  from  the  top  of  the  paper 
margin  to  begin  the  Y axis. 

- Increment  per  inch  on  the  X axis. 

- Increment  per  inch  on  the  Y axis. 


STARTX  - First  value  on  the  X axis. 

STARTY  - First  value  on  the  Y axis. 

LABELX  - X axis  label. 

LABELY  - Y axis  label. 

LEGEND  - Legend  for  the  curve. 


- File  where  the  X and  Y data  are  located. 


(7A10) 


ISA  10) 


(E10.3) 


(E10.3) 


(E10.3) 

(E10.3) 

(E10.3) 

(E10.3) 

(4A10) 

(4A10) 

(4A10) 


T 


NC0LX 

- Column  on  the  data  file  where 

the  X 

data  is 

(15) 

located. 

NC0LY 

- Column  on  the  data  file  where 

the  Y 

data  is 

(15) 

located. 

ISYM 

- Curve  plotting  symbol. 

(Rl) 

SC* LEX 

- Scale  factor  to  be  applied  to 

the  X 

data . 

(E10.3) 

SCALEY 

- Scale  factor  to  be  applied  to 

the  Y 

data. 

(E10.3) 

NOTE: 

1.  Cards  3 through  9 may  be  eliminated  by  using  option  12. 

2.  See  Paragraph  A.5.3.15  for  a sample  input. 

A. 5. 2. 15  Option  18  - Produce  a Cal  Comp  plot  of  the  time  history 

1. 


2. 


Input  the  following  cards  NPL0TS  times. 
5. 


1 


2-80 


1-5 

6-10 

11-80 

N0PT 

LEL0UT 

1-5 

6-10 

MUST 

NPRNT 

$PL0TD  NCURVS= , ANGLE3 , XPL0T= , YPL0T= , 


A-38 


(1 


3. 

1-10 

11-20 

-• 

K1 

NTAP 

Input  the  following  cards  if  NLIST  = 5. 

• # I 

n 

4. 

1-5 

6-10 

11-4 

li 

NPL0TS 

I 

LJ 

D 

D 

n 


LABEL 


where: 


N0PT 

- 

Option  number  (input  an  18  for  this  option). 

(15) 

LBL0UT 

- 

Option  description. 

(7A10) 

NLIST 

- 

File  where  the  plot  instructions  are  located. 

(15) 

NPRNT 

- 

Print  flag  to  display  input  (NPRNT=0,  no  print) 

(15) 

K1 

- 

Cal  Comp  plot  code. 

(A10) 

NTAP 

- 

File  where  Cal  Comp  writes  its  data. 

(A10) 

NPL0TS 

- 

Number  of  plots  to  be  produced. 

(15) 

I 

- 

Centimeter  plot  indicator. 

I = "C.M."  pen  movements  are  in  centimeter. 

(A4) 

NCURVS 

- 

Number  of  curves  per  plot. 

ANGLE 

- 

Plot  rotation  angle  in  degrees. 

XPL0T 

- 

Distance  in  inches  for  which  the  origin  of  the 
plot  is  to  be  moved  in  the  X direction. 

YPL0T 

- 

Distance  in  inches  for  which  the  origin  of  the 
plot  is  to  be  moved  in  the  Y direction. 

XEXIT 

- 

Final  pen  movement  in  the  X direction  in  inches. 

YEXIT 

- 

Final  pen  movement  in  the  Y direction  in  Inches. 

XLEGND 

- 

X coordinate  of  the  plot  legend  in  inches. 

YLEGND 

- 

Y coordinate  of  the  plot  legend  in  inches. 

X0RG 

- 

X value  at  which  axes  intercept 
(Default  = 0.0). 

Y0RG 

- 

Y value  at  which  axes  intercept 
(Default  = 0.0). 

LNGTHX 

- 

Length  of  X axis  in  inches. 

LNGTHY 

- 

Length  of  Y axis  in  inches. 

STARTX 


Initial  value  on  the  X axis. 


STARTY  - Initial  value  on  the  Y axis. 

DELX  - Increment  per  inch  on  the  X axis. 

DELY  - Increment,  per  inch  on  the  Y axis. 

IS0PT  - Bit  array  for  read  option.  See  Figure  A. 4 

BIT(O)  = 1 - Read  title  to  be  centered  over  plot. 

BIT(?)  = 1 - Pead  X axis  label . 

BIT(2)  - 1 - Read  Y axis  label . 

BIT(3)  = 1 - Read  legend  title. 

BIT(4)  = 1 - Read  coordinates  of  the  po’nt  and 
the  title. 

TITLE  - Four  card  plot  title. 

LABELX  - X axis  label. 

LABELY  - Y axis  label. 


LEGEND 


NDATA 

NCPI.X 


NC0LY 


- Legend  for  the  curve. 

- X coordinate  for  plot  title. 

- Y coordinate  for  plot  title. 

- File  where  the  plot  data  is  located. 

- Column  number  on  the  data  file  where  the  X 
data  is  located. 

- Column  number  on  the  data  file  w^ere  the  Y 
data  is  located. 

- Column  number  on  the  data  file  where  the  X 
elliptical  axis  data  is  located. 

- Column  number  on  the  data  file  where  the  Y 
elliptical  axis  data  is  located. 

- Column  number  on  the  data  file  where  the  cross 
correlation  data  is  located. 


(8A10) 

(4A10) 

(4A1.0) 

(4A10) 

(F10.0) 

(F10.0) 
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m 


LINTYP 


DASH 

NSYM 
SCALE X 

SCALEY 

SCALXX 

SCALYY 

LSTYLE 

LABEL 


Type  of  line  for  the  olot. 

LINTYP  = zero,  line  only  (no  symbol) 

LINTYP  = negative,  symbol  only  (no  lire) 

LINTYP  = positive,  both  line  and  symbol 

Length  of  plotting  segments  for  dashed  lines  in 
inches.  See  Table  A.i. 

Cal  Comp  symbol  code  (99  if  ellipse). 

Scale  factor  to  be  supplied  to  the  X data 
(Default  = 1.0). 

Scale  factor  to  be  applied  to  the  Y data 
(Default  = 1.0). 

Scale  factor  to  be  applied  to  the  X elliotical 
axis  (Default  = 1.0). 

Scale  factor  to  be  applied  to  the  Y elliptical 
axis  (Default  = 1.0) . 

Line  style.  See  Table  A.I.  (06) 

Legend  for  line  style.  (4A10) 


NOTE : 

1.  Cards  3 through  24  may  be  input  from  file  NLIST 
generated  by  option  13. 

2.  Ellipses  of  probability  may  be  plotted  bv  soecifying 
NSYM  = 99,  NCX,  NCY,  and  NCXY.  See  Figure  A. 3 

for  an  example  of  a plot  with  ellipses  of  probability. 

3.  See  Paragraph  A. 5. 3. 16  for  a sample  input. 

4.  See  Figure  A. 2 for  a diagram  of  the  plot  variables. 

A. 5.2. 16  Option  19  - Print  the  time  history  in  a blocked  form 


1. 

1-5 

6-10 

11-80 

N0PT 

LBL0UT 
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r 


! 


D 


1 


I 


U 

n 


i 


Input,  the  following  three  cards  as  many  times  as  needed  if  NLIST  = 5. 


1-5 

6-10 

11-15 

16-20 

21-70 

71-80 

I 

J 

K 

L 

LABEL 

SCALE 

Input  the  following  card  if  I < 0. 


8. 

1 

-40 

41-80 

FMTI 

FMT0 

9. 

1-10 

11-20 

21-30 

71-80 

R(I+J) 

R'I+J+1) 

R(I+J+2) 

• . . 

R(K+L) 

Input  the  following  two  cards  as  many  times  as  needed  if  NTAP  = 5. 


1-5 

6-10 

11-15 

16-20 

21-70 

71-80 

I 

J 

K 

L 

LABEL 

SCALE 

Input  the  following  card  if  I < 0. 


11. 

1-40 

FMT0 

where: 


N0PT 

- Option  number  (input  a 21  for  this  option) 

(15) 

LBL0UT 

- Option  description. 

(7A10) 

TITLE 

- Four  card  title  for  problem  definition. 

(8A10) 

NLIST 

- File  containing  the  R array  input  instructions 
for  the  integration  variables. 

(15) 

NTAP 

- File  containing  the  R array  output  instructions 
for  the  initial  integration  print. 

(15) 

NDATA 

- File  to  be  written  containing  the  time  history 
output.  Variables  set  up  by  option  7. 

(15) 

ii 
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Integration  option. 


NPRNT 


LABEL 


SCALE 


INT  = 1 - Fourth  order  Runge-Kutta. 

INT  = 2 - Third  order  Runge-Kutta. 

INT  = 3 - Second  order  Adams. 

- Print  the  time  history  as  a block  of  data. 

- Beginning  segment  number  of  the  R array. 

- Beginning  element  number  of  the  R array. 

- Ending  segment  number  of  the  R array. 

- Ending  element  number  of  the  R array. 

- Variable  names  for  the  elements  between 
I , J and  K,L. 

- Scale  factor  to  be  applied  to  R array  data 
between  I ,J  and  K,L. 

- Variable  input  format  for  inputting  data  into 
the  R array  for  card  9 above. 

- Variable  output  format  for  printing  data  from 
the  R array.  (Default  = 7F14.5). 


R array  elements. 


1.  Cards  7,  8,  and  9 are  input  through  Option  4.  The  branch 
to  Option  4 occurs  within  Option  21  and  is  invisible 
to  the  user.  Any  R array  element  may  be  input  through 
this  option.  To  terminate  the  read  for  cards  7,  8, 
and  9,  a blank  card  for  I,J,K,L  must  be  input. 


2.  If  NLIST  i 5,  cards  7 and  8 are  input  through  file 
NLIST.  Card  9 is  input  through  file  5 (INPUT). 


(4A10) 


(E10.3) 


(4A10) 


(4A10) 


(F10.3) 
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3.  Cards  10  and  11  are  input  through  option  5.  These  cards 
print  the  results  of  trim  and  the  call  to  Setup.  Any  R 
array  elements  may  be  printed  at  this  time.  To  terminate 
the  read  for  cards  10  and  11,  input  a blank  card  for 

I ,J,K,L. 

4.  If  NTAP  f 5,  input  cards  10  and  11  through  file  NTAP. 

5.  See  Paragraph  A. 5. 3. 18  for  a sample  input. 


A. 5.2. 18  Option  28  - Store  or  retrieve  R array  data  to  or  from 
a file  (TAPE28) 


1-5 

6-10 

11-80 

N0PT 

LBL0UT 

1-5 

K1 

1-5* 

6-10 

11-15 

16-20 

21-70 

71-80 

I 

J 

K 

L 

LABEL 

SCALE 

where: 

N0PT 

LBL0UT 

K1 


I 

J 

K 

L 


Option  number  (input  a 28  for  this  option). 
Option  description. 

Data  option  flag. 

K1  = 0 - Read  data  from  a file  (TAPE28). 

K1  / 0 - Write  data  to  a file  (TAPE28). 

Beginning  segment  number  of  the  R array  to  be 
read/written  from/ to  the  file. 

Beginning  element  number  of  the  R array  to  be 
read/written  from/ to  the  file. 

Ending  segment  number  of  the  R array  to  be 
read/written  from/ to  the  file. 

Ending  element  number  of  the  R array  to  be 
read/written  from/ to  the  file. 


(15) 

(7A10) 

(15) 


(15) 

(15) 

(15) 

(15) 
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LABEL  - Variable  names  appearing  between  segments  (5A10) 

I,J  to  K,L. 

SCALE  - Scale  factor  to  be  applied  to  the  R array.  (E10.3) 

NOTE: 

1.  Multiple  sets  of  card  3 may  be  input. 

2.  To  terminate  option  28,  a blank  card  for  I,J,K,L 
must  be  input. 

3.  See  Paragraph  A. 5. 3. 19  for  a sample  input. 


A . 5 . 3 Sample  Input  for  Each  Program  Option 

The  following  examples  are  sample  inputs  for  the  initialization  phase 
and  the  program  options. 

A. 5. 3.1  Sample  of  program  initialization 


SSEGMNT  NSEGI-7*  N$EGR«192#  NPRNT-1# 

NKl). 10*0, 0*19*19, 6* 

NR< l)-6, 17,17,26# 16.11 .37,66,9# 37, 

NR  ( 1?  )-19, If, 6*10, 0,6# 19 #0,0* 10# 

NR  ( 2 1 )■  1(5#  10#  10#10#10#0#6#6#6#6» 

NR  (31)  *6#  6#  6#  6#  6#  6#  6# 6*6*6# 
NR(91)-6#6»t,6»6»fc,6#6#0#289# 
NR(51)*269#289»99»16#50#50#75»75#75#  79# 
NR(61)»75#75*6#6»6#6#6#6#6#6# 

NR (71) »6 >6# 6# 6# 6# 6 #6* 6# 6# 6# 
NP(81)«36»29»92#99#28#1#1»1#1#29# 

NR (91 )-16#l#l#l*l#l«l#l»l#l» 

NR (101 >-l#l#l#l#l#l#12#6#6#6# 
NR)111)»6»6*6#6#6#6#6#6#6#6# 

NR  (121  )•  6*  6#  6*  6*  6*  6*  6*  6#  6#  6# 

NR (131) ■6*6*6*6*10# 10. 10# 10# 10*4# 

NR ( 191 ) *9#  9#  9# 


SENO 


T 

y 

u 

A. 5. 3. 

2 Sample 

of  Option  1 

^ 

□ 

1 

STORE 

R ARRAY  INSTRUCTIONS  ON  A FILE. 

n 

3 

6 

1 

6 

6 

WEIGHT. 

G.IX.IY. 12. IXZ 

u 

6 

7 

6 

9 

IN.  wn# 

LTOH 

7 

9 

7 

- 9 

TAUO 

• 

7 

28 

7 

29 

OHEGAV. 

KV 

7 

30 

7 

35 

TAUA.  Kl.  K2,  KY#  81.  82 

13 

1 

13 

6 

TAR  VI 

(V  TABLE  FOR  AERO  DATA  1 1.688 

1* 

1 

1* 

10 

TABV2 

(TINE  TABLE  FOR  OPEN  LOOP  INPUTS) 

20 

1 

20 

10 

TAPTC 

(OPEN  LOOP  THETA) 

120 

1 

120 

6 

TYV 

121 

1 

121 

6 

TYP 

122 

1 

122 

6 

TYR 

123 

1 

123 

6 

TYDA 

12* 

1 

12* 

6 

TYOR 

125 

1 

125 

6 

TLV 

126 

1 

126 

6 

TLP 

127 

1 

127 

6 

TLR 

128 

1 

128 

6 

TLOA 

129 

1 

129 

6 ' 

TLOR 

130 

1 

130 

6 

TNV 

IT 

131 

1 

131 

6 

TNP 

132 

1 

132 

6 

TNR 

• • 

133 

1 

133 

6 

TNDA 

is* 

1 

13* 

6 

TNOR 

135 

1 

135 

10 

TABUB 

136 

1 

136 

10 

TABWR 

137 

1 

137 

10 

TXAPP 

138 

1 

138 

10 

TZAPP 

139 

1 

139 

10 

TXSP 

1*0 

1 

1*0 

* 

TOAIN 

1*1 

1 

1*1 

* 

TDAOUT 

1*2 

1 

1*2 

* 

TDRIN 

1*3 

1 

1*3 

* 

TOROUT 

u 

BLANK  CARO  FOR  PE  AO  R 

n 

[J 

A. 5, 3. 

3 Sample  of  Option  2 

n 

2 

REAO 

INTO  I ARRAY 

I 

5 

5 

0 

1 

3 

1 

6 

NSTATES.  NPA.  NON.  NGN 

I 

11 

17 

* 

11 

1 

1 

7 

1 

10 

NPE.  NPX.  NOAW.  NOINIT 

[] 

1 

1 

0 

' 1 

BLANK 

CARO  FOR  READ  I 

i 

u 

A-49 

n 

— . — - — 

• ~ — - — - 

- .1 

A. 5. 3. 4 Sample  of  Option  3 


3 PRINT  DATA  FROM  THE  I ARRAY 

5 

1316  NSTATES#  NPA>  NON,  NGN 
1 7 1 10  NPE#  NP X#  NOAW#  NOINIT 

BLANK  CARD  TO  TERNINATE  OPTION  3 


A. 5. 3. 5 Sample  of  Option  4 


* 

5 

5 

READ  R 
0 

(CASE-PARTICULAR  INPU1 

8 

1* 

8 

1* 

VANS 

20. 

8 

15 

8 • 

15 

PS  IAMB 

150. 

6 

10 

8 

10 

V5  (SHIP  SPEED) 

20. 

8 

11 

8 

11 

PSI5  ((SHIP  HEADING) 

30. 

8 

7*. 71 

3 

8 

-7.012 

5 

UAS#  VAS#  WAS 

10. *95 

3 

-1. 

15 

3 

0. 

16 

PHI#  PSI 

5 

290. *1 

10 

5 

0. 

12 

SUHFYO#  SUMLO#  SUNNO 
0. 

2 

17 

2 

17 

YSP 

0. 

2 

0. 

1* 

? 

1* 

YAPP 

7 

10. 

36 

7 

2. 

37 

TL1#  TL2 

1* 

1. 

5 

16 

1. 

6 

GAINOA#  GAINDR 

A. 5. 3. 6 Sample  of  Option  5 


1.668 

.017*5329 

1.688 

.017*5329 

.017*5329 


5 

5 

2 1 

-8  1 

(6E20.10 ) 


PRINT  DATA  FRON  THE  R ARRAY 

2 5 UOOT#  WDOT#  QDOT#  TOOT#  XAPPOOT 

8 6 ALT#  VAIR#  UBAS,  VBAS#  WBAS#  ALPHA 

BLANK  CAPO  TO  TERMINATE  OPTION  5 
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A. 5. 3. 7 


Sample  of  Option  6 


6 

4 

f 

SETUP 

OIF 

FEREKTIAL  EOUATIONS 

■1 

2 

1 

11 

2 

17 

VDOT,PDOT*ROOT# YAPPDOT# PHIOOT# PSIDOT# YSPDOT 

12 

C 

15 

12 

18 

DA1D#DA2D#DP1D#DR2D 

BLANK  CARD  TO  TERMINATE  DERIVATIVES 

D 

3 

1 

11 

3 

17 

VB#  PB#  RB#  YAPP#  PHI#  PSI#  VSP 

12 

4 

11 

12 

15 

DA1#DA2#DR1#DR2 

BLANK  card  to  terminate  STATES 

C 

52 

Jk 

1 

1 

* 

51 

1 

PDOT,  P 

BLANK  CARD 

100 

l 

4 

99 

1 

PEO#  PE 

BLANK  CARO 

O 

105 

l 

103 

1 

PXHATD#  PXHAT 

BLANK  CARD#  END  OF  OPTION  6 

A. 5. 3,8  Sample  of  Option  7 


DESIGNATE  OUTPUT  VARIABLES 


1 

1 

1 

1 

TIME 

e 

1 

8 

2 

ALT#  VAIR 

3 

5 

3 

6 

XAPP#  ZAPP 

3 

15 

3 

15 

YAPP 

3 

17 

3 

17 

YSP 

3 

15 

3 

16 

PHI#  PSI 

12 

7 

1? 

8 

DA#  OR 

16 

3 

16 

5 

SIGOA#  SIGDR 

3 

11 

3 

11 

VB 

3 

1 

3 

2 

UB#  WB 

e 

27 

8 

29 

XM»  YM#  ZM 

8 

30 

e 

32 

PSIA,  thetaa#  PHIA 

8 

33 

e 

35 

SIGMAVX#  S I G** A VY#  SIGMAVZ 

6 

39 

8 

52 

UWINO#  UMEAN#  UR#  UAW 

6 

53 

8 

56 

VWIND*  VMEAN*  VR*  VAW 

8 

57 

8 

5C 

WVINO*  VMEAN#  WR*  WAV 

55 

1 

55 

17 

SIGMAS  OF  THE  STATES 

11 

9 

11 

9 

ERY 

16 

1 

16 

2 

SIGOAP#  SIGORP 

BLANK  CARD  TO  TERMINATE  OPTION  7 


A. 5. 3. 9 Sample  of  Option  10 


10 

■ 

CREATE 

A FILE 

COMPATIBLE 

FOR  OPTION 

IB 

* 

5 

S 0 

1.0 

2.0 

3.0 

5.0 

5.0 

6.0 

7.0 

0.0 

1.1 

2.1 

3.1 

5.1 

5.1 

6.1 

‘ 7.1 

B.l 

1.2 

2.2 

3.2 

5.2 

5.2 

6.2 

7.2 

a.  2 

1.3 

2.3 

3.3 

5.3 

5.3 

6*3 

7.3 

8.3 

1.5 

2.5 

3.5 

, 5.5 

5.5 

6*4 

7.5 

0.5 

1.5 

2.5 

3.5 

5.3 

5.3 

6.5 

7.5 

8.3 

A-51 


a 

A. 5. 3. 10  Sample  of  Option  11 


11 

2 

0 

STORE  TABBING  INSTRUCTIONS 

ii 

« 

10 

1 

TINE 

10 

2 

ALT  ' 

1 

10 

a 

VA-KTS 

.992A 

10 

8 

PHI 

97.3 

10 

9 

PSI 

97.3 

10 

A 

X APP 

n 

to 

6 

YAPP 

10 

7 

YSP 

• m 

10 

10 

99 

ERY 

n 

10 

1 

TINE 

10 

1 A 

VB 

10 

30 

WIND 

10 

31 

VNEAN 

10 

33 

VAN 

10 

2A 

SIGMAVY 

10 

10 

DA 

10 

11 

DR 

10 

12 

SIGOA 

lfi 

13 

SIGOR 

10 

10 

1 

TINE 

10 

38 

SI  GV 

u 

10 

39 

SIGP-OEG 

97.3 

10 

AO 

SIGR-DEG 

97.3 

10 

At 

SIGPMI 

97.3 

10 

A2 

SIGPSI 

97.3  ii 

10 

A3 

SIGY 

10 

20 

PSIA 

97.3 

10 

21 

THETAA 

97.3 

10 

22 

PHIA 

97.3 

9 

10 

1 

TINE 

10 

96 

SIGDAP 

10 

97 

SIGDRP 

• *. 

10 

12 

SIGDA 

10 

13 

SIGOR 

r t 

BLANK  CARO  TO  TERMINATE  OPTION  II 


D 

A-52 

1 


A. 5. 3. 11  Sample  of  Option  12 


U 


u 


12  STORE  PRINTER  PLOTTER  INSTRUCTIONS  ON  A FILE. 

3 

3 

VOLAP  CHECKOUT  - 16  FEB  78 
EFFECT  OF  AIRVAKE  ON  SIMPLE  MODEL 
UAH  VS  TIME 

1 0.0  7.0  3.0  1.0 

TIME  - SEC.  UAV 

UAH  IN  FEET  PER  SECONO.  10  1 40  U 

VOLAR  CHECKOUT  - 16  FEB  76 
EFFECT  OF  AIRWAKE  ON  SIMPLE  MOOEL 
HAW  VS  TIME 


1 0.0  0.0  3.0  1.0 

TIME  - SEC.  HAW 

HAW  IN  FEET  PER  SECONO.  10  1 *8  H 

VOLAR  CHECKOUT  - 16  FEB  78 
EFFECT  OF  AIRWAKE  ON  SIMPLE  MODEL 
VAW  VS  TIME 

1 0.0  7.0  3.0  1.0 

TIME  - SEC.  VAW 

VAW  IN  FEET  PER  SECOND.  10  14 A V 


FT/SEC. 

1.0  1.0 


FT/SEC. 

1.0  1.0 


FT/SEC. 

1.0  1.0 


A-53 


A. 5. 3. 12  Sample  of  Option  13 


13  STORE  CALCOMR  PLOTTER  INSTRUCTIONS  ON  A FILE. 

2 

CAL36  PLOT 

9 

SPLOTO  NCURVS-1#XPLOT-2.0#YPLOT-0.3,ANGIE-0.0»LNGTHX-5.0»STARTX-500.0# 
DELX— 100.0# X0RG-500. 0#LNGTHY«A.0# STARTY-O.O# DEL Y-2.0#YORG-0.0# 
ISOPT-22# SEND 

RANGE  FROM  SHIP  - FT. 

SIGMA  APP.  PATH  - FT. 

5.0  9.75 

VOLAR  CHECKOUT 
1*  FEB  78 
EFFECT  OF  AIRWAKE 
OS  SIMPLE  MODEL 

SCURVE  NDATA-10#NC0LX-21# NC PL Y-31#LINTYP-0» DASH -0.01# SC  ALEX -1.0# 
SCALEY-1.0#SEN0 

177777 

SPLOTO  NCURVS-2#XPL0T-0.0#YPL0T-5.5#XLEGN0-0.25#YLEGND"1.00# 

DEL Y- 5C.  0,1  SOP T-12# SEND 
ALTITUDE  - FT. 

5C  PERCENT  PROP.  ELLIPSE  X 5 
SCURVE  nC0LY-2#NSYM-99,LINTYP-10#NCX-23#HCY-2A#NCXY-53# 

SCALXX-5.0#  SC  AL YY-5 .0#  SENO 
177777  ACTUAL 

SCURVE  NCOLY-16,LINTYP-0#SCALEY--1.0#SENO 
170360  COMMANDEO 

SPLOTO  NCURVS-  l»XPL0T-7.5# YPLOT-0.0#DELY-2.0# ISOPT-A# SENO 
SIGMA  X A P - FT. 

SCURVE  NC0LY-23#SCALEY-1.0#SEN0 
177777 

SPLOTO  NCURVS-1 #XPLOT.O.C#YPLOT--5.5#OELY-2.0#SENO 
SIG**A  ZAP  - FT. 

SCURVE  NC0LY-24,SEN0 
177777 

SPLOTO  XPLOT-9.5#YPLOT-C.#OELY-2.  SENO 
SIGMA  THETA  - DEG. 

SC UP VE  NCOLY-22,  SCALEY-57.3  SENO 
177777 

S PL  CTO  XPLPT-C.0#YPL0T»5.5#STARTY--8.C#DELY-A.0#SEND 
THETA  - DEG. 

SCURVE  NCOLY-A  SENO 
177777 

SPLOTO  . XPlOT-7.5#YPL0T-0.#STARTY--.2#DELY-.l  SENO 
DE  - UNITS 

SCURVE  NCOLY-7#  SCALEY-1.  SENO 


SCURVE 

177777 

SPLOTO 


XPIPT-O.# YPL0T--5.5#STARTY-0.#0ELY-.05  SENO 
SIGMA  OE  - UNITS 
NCOLY-28  SENO 


SCURVE  NCOLY-28  SENO 

177777 

SPLOTO  XPL0T-9.5#YPL0T-C.#0ELY-10.0  SENO 
SIGMA  RPM  - RPM 
SCURVE  NCOLY-25  SENO 

177777 


A-54 


A. 5. 3. 13  Sample  of  Option  14 


14  TAB  THE  OUTPUT 

2 TAB  INSTRUCTIONS  STORED  ON  FILE  2 


A. 5. 3. 14  Sample  of  Option  15 


15  PUNCH  TIME  HISTORY 

3 


A. 5. 3. 15  Sample  of  Option  17 


17  PRODUCE  PRINTER  PLOTS 

5 1 

3 

VOLAR  CHECKOUT  - 16  FEB  76 


Effect 

OF  AIRVAKE  OK  SIMPLE 

MODEL 

* 

UAW  VS  TIME 

1 

0.0  7.0 

3.0  1.0 

TIME  - SEC. 

UAW 

FT /SEC 

UAW 

IN 

FEET  PER  SECOND. 

10  1 

4C  U 

1.0 

VOLAR  CHECKOUT  - 16  FEB 

78 

EFFECT 

OF  AIRVAKE  ON  SIMPLE 

MODEL 

WAW  VS  TIME 

1 

0.0  0.0 

3.0  1.0 

TIME  - SEC. 

WAW 

FT/SEC 

WAW 

IN 

FEET  PER  SECOND. 

10  1 

46  W 

1.0 

VOLAR  CHECKOUT  - 16  FEB 

78 

EFFECT 

OF  AIRWAKE  ON  SIMPLE 

MODEL 

VAW  VS  TIME 

1 

0.0  7.0 

3.0  1.0 

TIME  - SEC. 

VAW 

FT/SEC 

VAW 

IN 

FEET  PER  SECONO. 

10  1 

44  V 

1.0 

A-55 


A. 5. 3. 16  Sample  of  Option  18 


18  PRODUCE  CALCOMP  PLOTS 

5 1 

CAL3A  PLOT 

10 

SPLOTD  NCURVS-2»XPL0T-.75#YPL0T-0.0#ANGLE-0.0#LNGTHX-A.0»STARTX-0.0» 
delx-5.0# xprg-o.o# lngthy-a.o»starty-oo.o#dely-ao.o»yorg-o.o# 

ISOPT-fc  #XLEGND-2.0»YLEGND-1.0  *ENO 
TIME  - SEC. 

LATERAL  POS.  - M. 

SCURVE  DASH-.01,NDATA-10#NC0LX-1#NC0LY-6#LINTYP-0#SCALEX-1.0#SCALEY-.30A8  SEND 
177777 

SCURVE  DASH-.01#NDATA-10#NC0LY-7#SEND 
170360  COMMANDED 

SPLPTD  NCUR VS  - 1# X PLOT -0.0»Y PLOT- 5.0# l NGTHY-A.O# ST ARTY-0.0# DEL Y-  1.0#Y0RG-0.0» 
IS0PT-A,XLEGN0-2.0#YLEGND-A.5»S£ND 
SIGMA  LAT.  POS.  - h. 

SCURVE  DASH..01»NDATA«10»NC0LY-A3»SCALEX-1.0#SCALEY-.30A8  SEND 
177777 

SPLOTD  NCURVS-1#XP10T«5.3#YP10T--5.0#LNGTHY-A.0#STARTY— A.0»DELY-2.0# YORG-O.O# 
ISOPT-A,XLEGND«2.C#YLEGND-3.5  SEND 
BANK  ANGLE  - -DEG. 

SCUPVE  DASH-.01»NDATA-10#7CGLY-8»SCALEX-1.0#SCALEY-57.3  seno 
177777 

JPLPTD  NCURV$-1#XP10T-0.0#YPL0T-5.0#LNGTHY-A.0#STARTY-0.0#0ELY-  2.0# YORG-O.O# 
ISCPT-A,XLEGND-2.0#YLEGND«A.5#SEND 
SIGMA  PHI  - DEG. 

SCURVE  DASH-.01»NDATA-10#NC0LY-Al»SCALEX-1.0»SCALEY-57.3  SENO 
177777 

SPLOTD  NCURVS-1#XPL0T-5.5#YPL0T--5.0#LNGTHY-A.0# ST ARTY— 0. 2 # PEL Y-O.l# YORG-O.O# 
ISOPT-A#$END 

ROLL  CONTROL  - UNITS 

SCURVE  DASH-. 01# NO ATA-10*NC0LY-10*SCALEX-1,0»SCALEY-1,0» SEND 
177777 

S 0 LOTD  NCUR VS -1 » X PLCT-0.0#YPLnT- 5. 0#lNGTHY-A.0» ST  ARTY-0. 0# DEL Y-. 20# YORG-O.O# 
IS0PT-A,XLEGND-2.0#Yl EGND -A . 5# SEND 
SIGMA  ROLL  CTNTROL  - UNITS 

SCURVE  DASH-.01»NDATA-10.NCDLY-56,SCALEX-1.0»SCALEY-1.0,SEND 
177777 

SPLOTD  NCURVS-1#XPLOT-6.0»YPLOT--5.0#LNGTHY-A.O#STARTY-  -.02#DELY-.01# 
YOPG-O.O# ISPPT-A,SEND 
YAW  CONTROL  - UNITS 

SCURVE  DASH-. Cl»NDATA-10»NCnLY-ll#SCALEX-1.0»SCALEY-1.0# SEND 
177777 

SPLOTD  NCUR VS • 1. X PLOT- 0. 0# Y PLOT- 3.0# l NGTHY- A. 0# ST  ARTY-0.0# DEL Y-  .01# YORG-O.O# 
ISrPT-A#XLEGND-2.0»YLEGND-A.3#SEND 
SIGMA  YAW  CONTROL  - UNITS 

SCURVE  PASH-.0l#NDATA-10#NC0LY-57#SCALEX-1.0#SCALEY-1.0#SEN0 
177777 

SPLOTD  NCURV$-1.XPL0T-5.5#YPL0T--5.0»LNGTHY-A.0#STARTY-  -2.0#0ELY-1.0# 
YORG-O.O# ISOPT-A/SE NO 
LAT.  POS.  ERROR  - M. 

SCURVE  DASH-.01#NDATA-10#NC0LY-3 3# SC  ALEX-1.0# SC ALEY-. 30A8#  SENO 
177777 

SPLOTO  NC UR VS-1 » X PL OT -0. 0# YPl OT - 5. 0# L NGTHY -A, 0# ST ARTY-0. O/OELY-l.O# YORG-O.O# 
ISPPT-*»XLEGN0-2.0»YLEGND-A.5#SEND 
V AW  - M/SEC. 

SCURVE  DASH-.C1#NDATA-10#NCOLY-33#SCALEX-1.0#SCALEY-.30A8#SENO 
177777 


a 

fl 

.1 


a 


ri 


i 

i 


i 


a 

a 

a 


i 
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A. 5. 3. 17  Sample  of  Option  19 

19  PRINT  A FILE  IN  BLOCKED  FORK. 

» 

u 


A. 5. 3. 18  Sample  of  Option  21 


21  GENERATE  A TINE  HISTORY  BY  COVARIANCE  PROPAGATION. 

♦•♦♦•♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦♦♦♦♦A 


<M 

AV-BA 

** 

** 

35 

KT 

WOO, 

LAT/OIR  ** 

Li 

5 

5 

10 

3 

0 

1 

3 

1 

6 

OTIHE,  DTOUT,  0T0UT2,  THAX 

rr 

O.OI29 

0. 

2 

5.0  5.0 

BLANK  CARO  FOR  READ  R 

km 

e 

1 

e 

66 

PRINT  SEGMENT  « 

PLANK  CARD  FOR  WRITE  R 

A. 5. 3. 19  Sample  of  Option  28 


28 

STORE  R ARRAY  VALUES  ON 

TAPE20 

1 

51 

1 

FILE 

51 

TO  BE 
2764 

CREATEO 

P ARRAY 

99 

1 

99 

324 

PE  ARRAY 

103 

1 

103 

324 

PXHAT  ARRAY 

93 

1 

93 

9 

VU  ARRAY 

94 

1 

94 

36 

VY  ARRAY 

TERMINATE  OPTION  28 

BLANK  CARO  TO 

A. 5.4  Program  Output 

The  program's  output  is  divided  into  four  categories:  disk  storage, 
printed  output,  plotted  output  and  punched  output.  Before  the  time  history 
output  can  be  printed,  plotted,  or  punched,  it  has  to  be  stored  on  a file. 

A. 5.4.1  Disk  storage 

Disk  storage  is  divided  into  two  sections:  instructions  and  time  history. 
A. 5.4. 1.1  Instruction  storage 

Option  1 stores  I or  R array  instructions  on  a file.  Options  2 and  4 
have  the  capability  of  storing  I and  R array  instructions  with  the  data  on  a 
file  by  specifying  the  NPF  parameter.  Options  11,  12,  and  13  store  tabulating 
and  plotting  instructions  on  a file.  Option  28  stores  only  the  R array 
values  on  a file. 

A. 5. 4. 1.2  Time  history  storage 

Option  21  stores  the  time  history  output  (specified  by  Option  7}  on  a 
file.  Option  10  reads  input  and  creates  a file  of  the  same  structure  as  the 
time  history  file.  This  file  is  thus  compatible  with  the  program's  output 
features . 

A. 5. 4. 2 Printer  output 

The  time  history  output  can  be  printed  in  two  forms:  tabular  and  blocked. 
A. 5.4.2. 1 Tabular  output 

Once  the  time  history  has  been  generated  and  stored  on  a file  (A. 5. 4. 1.2), 
a tabular  print  of  the  output  can  be  produced  by  using  Options  11  and  14. 

Option  11  designates  which  elements  from  the  stored  file  (A. 5. 4. 1.2)  are  to  be 
printed  along  with  page  title  and  column  headings.  Option  14  produces  the 
tabular  print.  See  paragraphs  A. 5. 3. 10  and  A. 5. 3. 13  for  examples  of  Options  11 
and  14  respectively. 

A. 5. 4. 2. 2 Blocked  output 

The  contents  of  a file  can  be  printed  in  a blocked  form  by  using  Option 
19.  The  output  elements  are  not  titled  but  appear  in  the  same  order  as  they 
appear  on  the  file  being  dumped.  Option  21  has  the  capability  to  print  the 
time  history  in  a blocked  form  by  specifying  NPRNT  f 0.  This  is  to  be  used 
for  debugging  because  the  print  is  generated  every  DT0UT  seconds  until  job 


1 


A. 5. 4. 4 Punched  output 

After  the  time  history  has  been  placed  on  a file  (A. 5.4. 1.2),  punched 
cards  can  be  produced  by  using  Options  11  and  15.  Option  11  designates  which 
elements  from  the  stored  file  are  to  be  punched.  Option  15  produces  the 
punched  cards.  See  paragraphs  A. 5. 4. 10  and  A. 5. 4. 14  for  examples  of  Options  11 
and  i5  respectively. 

A. 6 Program  Flow 

An  overall  picture  of  program  flow  is  illustrated  in  Figure  A. 5. 

A. 7 Discussion  of  Control  Cards 

The  following  discussion  of  control  cards  is  sufficient  to  place  the 
VOLAR  program  online  and  execute  its  sample  problems. 

A. 7.1  Transferring  the  Data  from  Magnetic  Tape  to  Disk 

Table  A. 2 is  a description  of  how  the  tape  was  created. 


Table  A. 2 


TRACK 

LABEL 

CODE 

FORMAT 

PARITY 

DENSITY 

RECORD  TYPE 

BLOCKING 

NO.  FILES 

1 

2 

3 

4 

5 

6 


9 


None 

EBCDIC 

SI  (SCOPE  Internal) 
ODD  (Binary) 

1600  BPI 

S (Variable  length) 
Unblocked 


Classical  Pilot  Program  (UPDATE) 
Optimal  Pilot  Program  (UPDATE) 
Classical  Pilot  Data 
Optimal  Pilot  Data 
Classical  Pilot  Initialization 
Optimal  Pilot  Initialization 


Jl 

T 

i 

I 

ii 


A-60 


AIRWAKE  I DESCRIB  I TLU  I SOLVR 


The  following  card  setup  will  read  the  tape: 


J0B,NT1,CM15OOO,T1OO. 

ACC0UNT,cn,pw. 

FTN,A,B=lG0. 

DEFINE ,TAPE2=V0LAR1 ,TAPE3=V0LAR2 . 

DEFINE ,TAPE4=CASE 1 ,TAPE5=CASE2 . 

DEFINE, TAPE6=CASE1I ,TAPE7=CASE2I . 

REQUEST, TAPE 1,D=PE ,NT,P0=R,F=SI ,CV=EB,LB=KL ,VSN=5251 . 
LG0. 

7/8/9 

PR0GRAM  RDT AP E (TAPE1=0, TAPE2=0, TAPE 3=0, TAPE4=0, 
* TAPE5=O,TAPE6=O,TAPE7=O,0UTPUT) 

D I MENS I 0N  IBUF(30000) 

D0  100  J=2,7 

1 BUFFER  IN(1,1)  (IBUF(l) , IBUF( 30000 ) ) 

I F (UNIT ( 1 ) ) 5,100,5 

5 L=LENGTH(1 ) 

BUFFER  0UT(J , 1 ) ( I BUF ( 1 ) , I BUF ( L ) ) 

IN=UNIT (J ) 

GO  TO  1 

100  C0NTINUE 

ST0P 

END 

6/ 7/ 8/ 9 


A. 7. 2 Program  Execution 

The  following  card  setup  will  compile  the  Optimal  Pilot  program  and 
create  an  executable  file: 

J 0B , CM  1 50000 , T 500 . 

ACC0UNT,cn,pw. 

GET ,0LDPL=V0LAR2,TAPE28=CASE2I , DATA=CASE2 . 

UPDATE, F. 

DEFINE, LG0=V0LAR2E. 

FTN ,1 ,Sl ,R=3,A. 


A-62 


a 


II 

l) 


t J 


I 

1 

I 


tl 


a 

ii 

D 

1 


MAP, FULL. 

LOSET,PRESET=ZER0,LIB=CALC0MP . 

Li 

LG0.DATA. 

r i 

u 

7/8/9 

6/ 7/ 8/ 9 

A. 7. 3 Program  Modification 

The  following  card  setup  will  modify  one  of  the  routines  in  the  Optimal 
Pilot  program  and  execute  the  program. 

J0B.CM15OOOO.T5OO. 

ACCOUNT  ,cn,pw 

GET ,0LDPL=V0LAR2 ,0LD=V0LAR2E . 

UPDATE. 

FTN , I ,SL ,R=3,A. 

C0PYL. 

MAP, FULL. 

LDSET,PRESET=ZER0,LIB=CALC0MP. 

NEW. 

7/8/9 

UPDATE  CHANGE  CARDS  G0  HERE. 

7/8/9 

INPUT  DATA  CARDS  GO  HERE. 

6/ 7/ 8/ 9 


A. 8 Discussion  of  Input  Deck  for  the  Classical  Pilot  Model 

Simulation 

The  purpose  of  this  section  is  to  lead  the  prospective  user  of  VOLAR 
through  a sample  input  deck.  The  discussion  here  will  focus  on  the  necessary 
flow  of  options  to  generate  a time  history  and  point  out  any  special  input 
requirements  that  might  not  be  obvious  to  the  new  user.  This  input  deck  is 
identical  to  the  one  which  generated  the  AV-8A  classical  pilot  model  simula- 
tion and  plots  for  the  case  with  control  limiting.  After  studying  this 
section  the  user  should  run  this  case  on  his  computer. 

A schematic  of  the  job  setup  to  run  this  case  is  presented  in  Figure 
A. 6.  Figure  A. 7 is  a listing  of  the  input  deck.  The  circled  numbers 
provide  a cross  reference  to  this  test;  i.e.,  the  cards  following  © are 
discussed  below  in  Section  A. 8. 4. 

Throughout  the  discussion  of  this  input  deck, it  may  be  beneficial  to 
the  reader  to  refer  to  the  description  in  Section  A. 5 of  the  main  program 
and  its  options. 

A. 8.1  Namelist  Input 

The  first  step  is  to  provide  the  main  program  with  the  information 
required  to  bookkeep  the  I and  R arrays.  The  namelist  SEGMNT  establishes 
the  number  of  segments  in  each  array  and  the  number  of  elements  in  each 
segment.  The  format  for  NAMELIST  inputs  can  be  found  in  the  CDC  FORTRAN 
manual . 

The  namelist  SEGMNT  is  the  one  and  only  mandatory  input.  All  other 
inputs  are  either  user  selected  options  or  input  required  by  the  chosen 
option. 

For  this  case  the  I array  is  seven  segments  long,  and  the  R array  is 
composed  of  143  segments.  The  NPRNT  = 1 causes  the  I and  R array  segmenta- 
tion information  to  be  output.  The  lengths  of  the  segments  are  consistent 
with  the  I and  R arrays  presented  in  Appendix  E. 

A. 8. 2 Option  2 

Integer  values  that  define  system  size  and  control  the  selection  of 
certain  calculations  are  input  here.  These  parameters  are  defined  in 
Appendix  E. 

Although  NPE  and  NPX  are  dimensions  of  matrices  used  only  in  the  optimal 
pilot  case, they  are  given  a value  of  one  (1)  here.  The  reason  is  the 
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integration  subroutines  expect  to  integrate  a set  of  means  and  three  matrices. 
By  giving  NPE  and  NPX  a small  value,  a negligible  increase  in  run  time  is 
traded  fcr  the  simplicity  of  having  only  one  set  of  integration  subroutines. 
This  keeps  the  number  of  problem  dependent  subroutines  to  an  absolute 
minimum. 

This  run  will  include  the  ship  airwake  (NOAW  = 0)  and  is  not  an 
initialization  run  (NOIN IT  = 1). 

A. 8. 3 Option  4 

Required  values  of  the  R array  are  input  next.  At  Vought  fixed  geometric 
and  aerodynamic  properties  are  input  here.  Case-dependent  information  is 
input  later  in  the  deck.  At  the  beginning  of  job  execution,  the  'LDSET, 

PRESET  = ZERO'  control  card  initializes  core  to  zero.  Therefore,  parameters 
not  initialized  through  Option  4 begin  with  a value  of  zero. 

Note  that  the  inertias  are  input  here  as  1 (obviously  not  the  true  value). 
However,  they  are  not  required  to  generate  accelerations  (Vought  used 
dimensional  derivatives)  but  inertias  do  appear  in  SETUP  in  a denominator. 

This  is  all  related  to  the  calculation  of  forces  and  moments  in  DEQU,  where 
the  total  forces  and  moments  that  drive  the  mean  equations  must  be  obtained 
from  a reference  condition  and  the  perturbation  forces  and  moments. 

Tables  defining  the  independent  variable  arrays,  open  loop  inputs, 
aerodynamic  data  and  nonlinear  elements  are  also  input  here. 

The  weight  and  aerodynamic  derivatives  are  from  Reference  15.  The  ship 
geometry  inputs  represent  a DD-963.  Feedback  gains  and  time  constants  are 
for  the  classical  pilot  model  discussed  in  Section  6.  The  values  of 
TABUB,  TABWB,  TXAPP,  TZAPP,  and  TXSP  were  derived  from  the  task  description. 
These  tables  provide  the  reference  longitudinal  values  required  by  the 
lateral/directional  equations.  The  final  four  tables  define  the  nonlinear 
functions  (control  limiters)  discussed  in  Section  4. 

A. 8. 4 Option  6 

Here  the  user  specifies  the  variables  that  are  to  be  integrated  and  the 
variables  that  contain  the  integrated  values.  This  option  allows  one  to  use 
any  variable  in  blank  COMMON  as  an  integrand  or  integrated  value. 


Note  that  the  setup  includes  the  matrices  RED,  PE,  PXHATD,  and  PXHAT. 
These  are  associated  with  the  optimal  pilot  model  and  not  required  by  the 
classical  model.  In  Section  A. 8. 2 it  was  noted  that  these  are  specified  as 
1 by  1 matrices  (NPE  and  NPX)  and  included  in  this  setup  to  avoid  the  need 
for  two  sets  of  integration  subroutines. 

A. 8. 5 Option  7 

The  variables  specified  here  appear  on  the  output  file.  Their  selection 
is  totally  up  to  the  user.  All  of  the  output  features  (tabbing,  plotting, 
etc.)  access  the  file  containing  the  Option  7 specified  variables. 

A. 8. 6 Option  11 

The  variables  to  be  tabbed  are  user  defined.  Here  it  was  attempted  to 
tabulate  all  the  key  parameters  for  interpreting  the  performance  of  a given 
task. 

A. 8. 7 Option  4 

The  case-particular  input  goes  here.  The  order  is  not  important. 

The  first  two  terms  in  this  input  are  ambient  wind  magnitude  and  direc- 
tion. The  ambient  wind  gets  into  the  calculation  of  airspeed  and  wind  over 
deck. 

The  ship's  speed  and  heading  are  also  input.  The  ambient  wind  and  ship 
conditions  were  selected  to  give  35  knots  wind  over  deck  comming  at  30  degrees 
to  port.  This  value  is  a 'breakpoint'  in  the  ship  airwake  model  data  tables. 

The  initial  body  axis  components  of  airspeed  are  necessary  to  initialize 
the  airplane  inertial  speed.  Aside  from  this  exception,  airspeeds  are  calcu- 
lated in  DEQU  from  inertial  speeds.  This  particular  representation  was  chosen 
because  it  was  felt  that  most  often  the  airplane  flight  conditions  would 
be  specified  by  stating  an  airspeed,  angle-of-attack,  and  sideslip.  If  this 
is  undesirable,  the  user  can  modify  DEQU  to  meet  his  needs.  The  mean  com- 
ponents of  the  airwake,  that  enter  into  the  calculation  of  airspeed,  may  be 
obtained  by  requesting  a zero  second  time  history  after  specifying  the  air- 
plane and  ship  positions  in  space. 

The  initial  aircraft  bank  angle  is  -1  degree,  and  the  aircraft's  heading 
is  zero  degrees.  The  bank  angle  is  required  to  balance  out  the  initial  side 
force  on  the  airplane.  The  initial  side  force  is  due  to  the  non-zero  value 
of  sideslip  (VAS  + 0). 


The  value  of  SUMFYO  was  obtained  from  Yy.  The  initial  values  of  SUMLO 
and  SUMNO  are  zero.  It  is  assumed  that  the  roll  and  yaw  controls  would  trim 
out  the  moments  due  to  and  L^  so  that  the  total  moments  on  the  airplane 
are  zero.  SUMFYO,  SUMLO,  and  SUMNO  are  total  force  and  moments  on  the  air- 
plane and  are  required  in  subroutine  SETUP  to  initialize  YO,  LO,  and  NO. 

The  initial  displacements  of  the  ship  and  airplane  along  the  y earth 
axis  are  input  as  zero.  The  x-axis  displacements  of  the  ship  and  airplane, 
and  the  z axis  displacement  of  the  airplane  were  previously  input  in  table 
form. 

A. 8. 8 Option  28 

This  option  reads  data  stored  from  a previous  run.  It  is  used  here  to 
initialize  the  P matrix. 

In  the  previous  run,  NOIN  IT  was  set  to  -1  (all  other  input  discussed 
to  this  point  remains  unchanged).  This  froze  the  means  and  integrated  only 
the  covariance  matrix.  When  NOIN  IT  = -1,  Option  28  should  follow  Option  21. 
Refer  to  Section  A. 5 for  using  Option  28  in  a BUFFER  IN  or  BUFFER  OUT  mode. 

This  initialization  is  necessary  to  start  the  P matrix  at  realistic 
values.  Because  the  airplane  did  not  suddenly  appear  at  a point  16  seconds 
away  from  the  ship,  beginning  with  P = 0 is  unrealistic.  The  airplane  flew 
through  the  atmosphere  (and  its  disturbances)  to  get  to  the  point  in  the 
trajectory  where  the  VOLAR  simulation  picks  it  up.  The  question  then  becomes 
one  of  determining  what  is  a realistic  method  to  use  for  initializing  the 
covariance  matrix. 

Subroutine  SOLVR  provides  a method  for  initializing  P.  SOLVR  calculates 
the  steady-state  (i.e.  t -*•  ■»)  value  of  P.  This  solution  would  provide  an 
upperbound  on  the  uncertainties  (i.e.  the  value  of  P).  However,  SOLVR  will 
not  work  for  the  optimal  pilot  cases  because  there  are  three  inter-related 
covariance  matrices,  and  it  is  programmed  to  solve  only  one.  Also,  one  might 
argue  that  t ~ statistics  may  not  be  realistic  either.  Note  that  if  the 
system  is  unstable,  P -*•  ® as  t -*•  ». 

The  approach  selected  here  was  to  integrate  the  P matrix  for  some  period 
of  time.  Five  seconds  was  selected.  For  a stable  system,  the  P matrix  will 
approach  its  steady  state  value  fairly  rapidly.  This  is  something  the  user 
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must  evaluate  on  a case  by  case  basis.  For  the  runs  presented  here, 

5 seconds  was  quite  adequate.  One  more  comment  is  in  order.  When  making 
initialization  runs.it  is  often  necessary  to  use  a lower  integration  incre- 
ment and/or  a higher  order  integration  scheme  than  in  the  simulation  runs. 

The  reason  is  the  elements  of  the  P matrix  often  undergo  wild  perturbations, 
before  settling  out,  when  all  the  elements  are  initialized  to  zero. 

A. 8. 9 Option  21 

This  option  generates  time  histories  of  the  means  and  covariances.  The 
example  presented  here  asks  for  a 20  second  time  history  using  Adams  integra- 
tion. The  covariance  matrix  will  be  printed  out  at  the  initial  time  point 
(t  = 0)  and  every  20  seconds.  The  output  file  is  TAPE10. 

After  the  call  to  SETUP,  before  the  integration  process  begins, 
the  elements  of  segment  eight  will  be  output.  This  was  done  to  check  on  the 
parameters  associated  with  the  airwake  calculations.  It  is  not  necessary  to 
output  anything  here.  But  often,  the  user  desires  to  know  some  parameter 
values  but  does  not  want  to  take  up  space  on  the  output  file.  This  is  a 
convenient  place  to  look  at  the  values  of  such  parameters. 

For  a given  stepsize,  Adams  integration  is  three  times  faster  than  the 
third  order  Runge-Kutta  and  four  times  faster  than  the  fourth  order.  However, 
the  Adams  integration  stepsize  will  normally  have  to  be  lower  than  the  Runge- 
Kutta  stepsizes  to  obtain  the  same  degree  of  accuracy.  We  have  found  the 
stepsize  change  does  not  offset  the  speed  advantage;  and  for  production 
running,  Adams  integration  is  preferable.  As  the  user  modifies  these  routines 
and/or  treats  new  problems,  he  should  spot  check  his  integration  to  confirm 
the  accuracy  of  the  simulation.  Integration  problems  normally  manifest 
themselves  as  time  history  divergence,  indefinite  parameter  values,  or 
infinite  parameter  values. 

A. 8. 10  Option  14 

This  option  causes  the  output  file  to  be  tabbed  according  to  the 
instructions  stored  by  Option  11. 
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A. 8. 11  Option  18 

For  some  runs,  the  user  may  desire  Cal  Comp  plots  as  well  as  tabbed 
output.  The  final  option  in  this  setup  directs  the  plotting  of  15  parameters 
from  the  output  file.  Normally  one  does  not  plot  in  the  same  run  that 
generates  the  output  file.  It  is  difficult  to  know  a priori  good  scales  for 
the  plots.  The  scales  presenting  these  data  were  selected  after  having 
seen  a tab  output  of  the  desired  parameters. 

A good  description  of  the  input  required  by  the  plot  package  has 
already  been  presented  in  Section  A. 5.2.  The  best  approach  to  setting  up  plot 
instructions  is  to  lay  out  the  desired  axes  on  regular  graph  paper,  and 
use  the  grids  to  develop  the  distances  for  pen  movement,  etc. 

A. 9 Discussion  of  an  Input  Deck  for  the  Optimal  Pilot  Model 

Simulation 

The  purpose  of  this  section  is  to  lead  the  prospective  user  of  VOLAR 
through  a sample  input  deck  for  the  optimal  pilot  model  simulation.  The  user 
should  be  familiar  with  Section  A. 5 which  discusses  the  main  program  and  its 
options.  The  discussion  here  will  focus  on  the  necessary  flow  of  options  to 
generate  a time  history  and  point  out  any  special  inputs  that  might  not  be 
obvious  to  the  new  user.  This  input  deck  is  identical  to  the  one  which 
generated  the  20  second  sea  state  5 time  histories  appearing  in  Figure  26. 

A schematic  of  the  job  setup  to  run  this  case  is  presented  in  Figure  A. 8. 
Figure  A. 9 is  a listing  of  the  input  deck.  The  circled  numbers  were  added 
to  provide  a cross  reference  to  this  text;  i.e.  the  cards  following  @ are 
discussed  below  in  Section  A. 9. 4. 

After  studying  this  example,  the  user  should  run  this  case  on  his 
computer. 

A. 9.1  Namelist  Input 

The  first  step  is  to  provide  the  main  program  with  the  information 
required  to  bookkeep  the  I and  R arrays.  The  namelist  SEGMNT  establishes 
the  number  of  segments  in  each  array  and  the  number  of  elements  in  each 
segment. 

The  namelist  SEGMNT  is  the  one  and  only  mandatory  input.  All  other 
inputs  are  either  user  selected  options  or  input  required  by  the  chosen 
option. 
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5SEGMNT  NSEGI-7,  NSEGR-143,  NPRNT»1* 

Ni<n*io, 0,0,14,14, 6, 

NRI1I.6, 17,17,26. 18# 11, 37,66,9, 37, 
NR<11)»14»18*6,10»0»6»19,0,0, 10, 
NR<21)«10,10,10,10, 10,0,6,6,6.6. 
NR<31)»6»6,6»6*6»6*6»6»6»6, 

NR < 41 )«6» 6, 6, 6, 6, 6, 6, 6, 0,289, 

NR < SI  I -2 89, 289, 66, 16, 50, 50,75, 75, 75, 75, 
NR  <61 ) >75, 75, 6, 6, 6, 6, 6, 6, 6,  6, 
NR(71I»6,6,6, 6, 6, 6, 6, 6,6,6, 

NR <81 1 Oft, 24, 42, 49, 28, 1,1, 1,1, 24, 

NR < 91 )>16, 1,1,1, 1,1, 1,1, 1,1, 

NR < 101  l«l, 1,1. 1,1, 1,12, 6, 6, 6, 
NR<111)*6*6»6*6»6»6»6»6,6*6» 

NR  < 121 )a6,6,6,6,6,6, 6, 6, 6, 6, 

NR <131 >-6,6, 6, 6, 10, 10, 10, 10, 10. 4, 

NR  < 141 1*4, 4, 4, 
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Here  the  I array  is  seven  segments  long  and  the  R array  is  composed  of 
119  segments.  The  R array  is  shorter  than  the  previous  example  because  at 
Vought  we  actually  ran  the  longitudinal  example  first  then  later  expanded 
the  R array  to  include  the  additional  lateral/directional  terms.  The  new 
parameters  were  mostly  added  to  the  end  of  current  segments  or  new  segments 
generated  for  them.  The  expanded  R array  will  work  for  this  case  provided 
the  segment  length  information  is  adjusted  to  the  correct  values  here.  The 
lengths  which  need  to  change  are  those  which  define  matrix  sizes.  These 
have  the  word  'varies'  under  the  column  'element'  in  Section  E.2.  Note  that 
in  the  previous  listing,  segment  51  (the  P matrix)  was  289  long  and  here  it 
is  1521  long.  This  is  because  P went  from  a 17  by  17  matrix  in  the  classical 
pilot  example  to  a 39  by  39  matrix  in  the  optimal  pilot  example. 

One  should  compare  blank  COMMON  for  the  two  listings  in  Appendices  B 
and  C.  Note  that  they  are  essentially  the  same  except  for  the  new  segments 
and  the  added  terms  on  the  end  of  segments;  then,  verify  the  segment  lengths. 

Remember  that  the  main  program  uses  the  information  input  by  SEGMNT  to 
establish  a one  to  one  relationship  with  the  elements  of  blank  COMMON.  If 
the  segment  lengths  are  incorrectly  specified,  the  program  will  self  destruct; 
or  at  best,  not  run. 


A. 9. 2 Option  2 

Integer  values  that  define  system  size  and  control  the  selection  of 
certain  calculations  are  input  here.  The  elements  of  the  I array  are  defined 
in  Appendix  E. 

There  are  11  mean  equations  to  be  integrated.  The  three  covariance 
matrices,  P,  PE,  and  PXHAT  are  n by  n matrices  where  n is  39,  18,  and  18 
respecti vely . QM  is  a 15  by  15  matrix  (NQM  = 15)  and  GM  is  32  (NGM  = 32) 
by  15  (NQM).  This  run  will  include  the  ship  airwake  calculations  (N0AW  = 0) 
and  is  not  an  initialization  run  (N0INIT  = 1). 


A. 9. 3 Option  4 

Required  values  of  the  R array  are  input  next.  At  Vought,  we  normally 
input  fixed  geometric  and  aerodynamic  properties  here;  inputting  case  dependent 
information  later  in  the  deck.  At  the  beginning  of  job  execution,  the  control 
card  ' LDSET , PRESET  = ZER0.'  initializes  core  to  zero;  therefore  parameters 
not  initialized  through  an  Option  4 begin  with  a value  of  zero. 
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Note  that  the  value  of  IY  is  1;  obviously  not  a realistic  inertia.  The 
1 is  used  because  the  inertias  are  not  needed  to  generate  accelerations 
(because  of  dimensional  derivatives)  but  inertia  does  appear  as  a denominator 
in  SETUP.  The  numerator  term  is  zero  so  anything  in  the  denominator  will 
work.  This  is  all  related  to  the  calculation  of  forces  and  moments  in  DEQU, 
where  the  total  forces  and  moments  that  drive  the  mean  equations  must  be 
obtained  from  a reference  condition  and  perturbation  forces  and  moments. 

All  the  tables  that  define  the  open  loop  inputs,  dimensional  derivatives, 
and  optimal  pilot  feedback  gains  are  input  here.  They  are  consistent  with 
the  independent  variable  arrays  TABV1  and  TABV2. 

The  ship  geometry  terms  are  for  a DD-963  and  the  aerodynamic  derivatives 
are  from  Reference  15.  The  optimal  pilot  feedback  gains  were  determined 
consistent  with  the  derivation  in  Section  7. 

Note  the  breakpoints  in  TABV1  (segment  13,  elements  1-6).  The  only 
two  'legitimate'  points  are  35.  and  56.6  knots.  The  data  tables  that  are  a 
function  of  airspeed  are  input  with  the  same  value  at  0 as  35  knots  and  the 
same  values  for  60.,  61.,  and  62.  knots  as  56.6  knots.  This  was  done  to 
fill  up  the  tables,  and  in  the  case  of  0 knots,  to  keep  the  computer  from 
extrapolating  to  speeds  lower  than  35  knots.  Here  the  simulation  uses 
35  knots  data  for  any  airspeed  that  drops  below  35  knots.  Thirty-five  knots 
is  the  lowest  airspeed  the  airplane  should  have  (35  is  WOO  speed)  and  the 
56.6  knots  was  a value  initially  selected  for  the  beginning  of  the  time 
history.  The  pilot  feedback  gains  and  the  dimensional  derivatives  are 
assumed  to  vary  linearly  between  these  speeds.  To  cut  down  on  run  time,  the 

■ 

airplane  was  moved  closer  to  the  ship,  to  start  the  simulation,  and  the 

i m 

initial  airspeed  became  44.7  knots.  There  was  no  requirement  to  change  the 
data  tables  because  the  table  look-up  routine  will  give  values  for  44.7  knots 
just  as  it  will  for  all  airspeeds  between  35  knots  and  56.6  knots. 

A. 9. 4 Option  6 

The  user  next  specifies  the  variables  that  are  to  be  integrated  and  the 
locations  which  contain  the  integrated  values.  There  is  no  restriction  on  the 
number  of  variables  one  desires  to  integrate  as  long  as  the  fourth  and  fifth 
I array  segments  are  long  enough  to  keep  all  the  R array  indices. 

0 
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Also  note  that  the  number  of  mean  equations  should  equal  NSTATES  (I  array 
segment  1,  element  3). The  integration  subroutines  will  only  integrate  the 
first  NSTATES  derivatives  specified  in  the  first  part  of  Option  6.  When  setting 
up  matrices  to  be  integrated  it  is  only  necessary  to  define  the  R array 
location  of  the  first  element.  The  number  of  -Matrix  elements  integrated  is  deter 
mined  internally  from  the  matrix  size  information  read  in  previously  (Option  2). 

Option  6 expects  location  information  for  all  three  pairs  of  matrices 
(PDOT  and  P,  etc.)  whether  the  user  intends  to  integrate  these  or  not.  A 
good  example  of  this  is  the  classical  pilot  model.  Refer  to  Sections  A. 8. 2 
and  A. 8. 4 for  additional  comments  on  this  matter. 


A. 9. 5 Option  7 

The  variables  specified  here  appear  on  the  output  file.  Their  selection 
is  totally  up  to  the  user.  All  of  the  output  features  (tabbing,  plotting, 
etc.)  access  the  file  containing  the  Option  7 specified  variables. 

A . 9 . 6 Option  11 

The  variables  to  be  tabbed  are  user  defined.  Here  we  have  attempted  to 
tabulate  all  the  key  parameters  for  interpreting  the  performance  of  a given 
task. 


A. 9. 7 Option  4 

The  case-particular  input  goes  here;  i.e.,  the  things  necessary  to 
define  the  environment  (wind  speed,  etc.)  and  the  initial  conditions  of  the 
ai rplane  and  ship. 

The  initial  components  of  airspeed  (segment  8,  elements  3-5)  and 
initial  forces  and  moment  (segment  5,  elements  1-3)  are  required.  The 
airspeed  requirement  could  be  relaxed,  if  desired,  according  to  the  discussion 
of  Section  A. 8. 7.  If  the  equations  describing  the  means  in  DEQU  were 
perturbation  equations  or  it  was  not  desired  to  integrate  the  means  then  the 
initial  forces  and  moment  requirement  would  not  longer  be  valid. 

The  ship  motion  inputs  KZS,  ZETAZS,  and  WNZS  were  obtained  from  Section  9 
of  this  report.  The  value  of  KZS  is  the  value  of  K^  from  the  table  multiplied 
by  wn  squared.  ZETAZS  and  WNZS  are  the  damping,  c»  and  natural  frequency, 
o)n,  presented  in  Table  9 of  Section  9. 


A. 9. 8 Option  28 

This  option  reads  data  stored  from  a previous  run.  It  is  used  here  to 
initialize  the  matrices  P,  PE,  PXHAT,  VU,  and  VY.  A good  discussion  of  the 
reason  for  this  initialization  appears  in  Section  A. 8. 8 and  will  not  be 
repeated  here. 

The  important  thing  is  that  this  setup  expects  a TAPE28  to  be  ATTACH 'ed 
(GET'ed)  to  this  job;  and  for  this  file  to  contain  numerical  values  for  all  the 
elements  of  the  above  mentioned  matrices. 

To  generate  the  TAPE28  file  for  initialization  the  user  makes  the 
following  changes  to  the  run  deck.  In  Option  2 the  value  of  NSTATES  is 
reduced  by  one  and  the  value  of  N0 IN  IT  is  set  to  -1.  The  cards  setting 
RCD  as  an  integrand  and  RC  as  the  integrator's  output  should  be  removed 
from  Option  6.  In  one  of  the  Option  4's,  input  RCD  (segment  11,  element  2) 
with  the  value  -37.28  (its  t = 0 value).  Run  a 5 second  time  history  using 
third  or  fourth  order  Runge-Kutta  integration.  An  integration  stepsize  of 
.025  is  adequate.  Following  the  time  history  (Option  21)  call  Option  28  and 
store  the  desired  matrices;  remembering  to  include  the  necessary  control 
cards  to  catalog  this  file  for  future  use. 

A. 9. 9 Option  21 

This  option  generates  time  histories  of  the  means  and  covariances.  The 
example  presented  here  asks  for  a 20  second  time  history  using  Adams  integra- 
tion. The  output  of  MISCAL  will  be  printed  at  the  initial  time  (t  = 0)  and 
every  20  seconds.  The  output  file  is  TAPE10. 

After  the  initial  call  to  setup,  before  the  integration  process  begins 
the  elements  of  segment  eight  will  be  output.  This  was  done  to  check  on  the 
parameters  associated  with  the  airwake  calculations.  It  is  not  necessary 
to  output  anything  here.  But  often,  the  user  desires  to  know  some  parameter 
values  but  does  not  want  to  take  up  space  on  the  output  file.  This  is  a 
convenient  place  to  look  at  the  values  of  such  parameters. 

A. 9. 10  Option  14 

This  option  causes  the  output  file  to  be  tabbed  according  to  the 
instructions  stored  by  Option  11. 
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A . 10  A Discussion  of  the  Programmed  Equations  in  Subroutines  DEQU  and  PROP. 

DEQU  and  PROP  are  the  primary  subroutines  whereby  the  user  defines  his 
problem  for  analysis.  Their  importance  warrants  a more  detailed  explanation 
than  previously  Dresented  in  Section  A. 3. 

This  discussion  addresses  the  "classical  pilot  model"  example.  Once  the 
user  is  familiar  with  this  example,  transition  to  the  more  complex  "optimal 
pilot  model"  example  will  be  simple.  Sub-section  titles  appearing  below  in 
quotation  marks  are  the  same  "Titles"  preceding  a group  of  calculations  in 
the  subroutines.  This  provides  the  user  with  a one-to-one  correspondence 
between  this  section  and  the  listings  of  Appendix  B. 

DEQU 

Subroutine  DEQU  provides  the  main  program  with  the  differential  equations 
for  the  mean  values  of  all  state  and  control  variables.  It  also  provides  the 
zero  subscripted  (reference)  values  of  these  variables  which  define  the 
conditions  about  which  the  motions  are  linearized. 

"Look-Up  Values  of  Required  Longitudinal  Parameters" 

Recall  that  the  classical  pilot  model  example  (Volume  I,  p.  41)  employs 
lateral/directional  airplane  dynamics.  Nevertheless,  the  associated  equations 
require  some  longitudinal  parameters.  As  an  example,  consider  the  equation 
for  rate  of  change  of  roll  attitude: 

♦ = P + Q sin(<J>)  tan(0)  + R cos  (4>)  tan(0) 

Q and  0 are  longitudinal  parameters.  Required  to  generate  <J>. 

This  section  of  coding  provides  the  necessary  longitudinal  parameters  as 
a function  of  time.  Values  of  the  longitudinal  parameters  for  the  reference 
trajectory  are  stored  in  the  appropriate  tables;  at  the  desired  time  instant, 
values  are  determined  by  table-look-up.  When  performing  an  initialization  run 
the  table- lock-ups  are  made  only  once. 

"Miscellaneous  Calculations" 

Calculations  performed  here  include: 

(a)  0 = 0C.  This  sets  airplane  pitch  attitude  equal  to  the  pitch 
attitude  froi..  the  table-look-up  of  the  previous  section. 


CJ 
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(b)  Ship  vertical  position  is  set  to  zero  (ZSP  = 0). 

(c)  Perturbation  values  for  roll  rate  and  yaw  rate  are  evaluated;  for 
example,  p = P - PQ,  where  p is  the  perturbation,  P is  the  total 
roll  rate  and  P0  is  the  reference  value.  The  reference  values  are 
updated  at  the  end  of  each  integration  increment  as  part  of  the 
quasi  linearization  process. 

(d)  The  remaining  calculations  are  performed  to  minimize  the  number  of 
calls  to  the  trigonometric  functions  used  throughout  the  remainder 
of  the  subroutine. 

"Altitude" 

Due  to  the  standard  axis  convention,  altitude  is  "minus  z". 

"Airplane  With-Respect-To  Ship,  Earth  Axis" 

Both  the  ship  and  the  airplane  are  referenced  to  a fixed  earth  axis 
system.  These  equations  calculate  the  x,  y,  and  z components  distance  between 
the  aircraft  c.g.  and  the  ship  c.g.  in  that  earth  axis  system. 

"Airplane  Wi th-Respect-To  Ship,  Ship  Wind  Axis" 

The  airwake  model  requires  the  coordinates  of  the  airplane  relative  to 
the  ship  to  be  in  ship  wind  axes.  These  equations  transform  the  earth  axis 
distances  to  ship  wind  axes.  The  parameters  DXM145,  etc.  were  calculated  in 
subroutine  SETUP. 

"Euler  Angles  for  Airwake  Model" 

These  are  Euler  angles  of  the  airplane  with  respect  to  the  ship  wind 
axes  and  are  necessary  inputs  to  the  airwake  model  eauations.  They  are 
simplified  greatly  because  the  airwake  model  used  does  not  pitch  and  roll  with 
the  ship. 

"Ship  Airwake" 

If  ship  airwake  is  desired  (NOAW  = 0)  subroutine  AIRWAKE  is  called. 
Following  the  call  to  AIRWAKE,  the  airwake  related  parameters  are  available 
in  COMMON . 

"Airplane  Body  Axis  Components  of  Ambient  Wind" 

These  equations  transform  the  ambient  wind  velocity,  which  is  given  in 
earth  axes,  to  the  aircraft  body  axes.  Wind  velocities  are  defined  consistent 
with  the  earth  axis  reference,  i.e.  a positive  x-axis  wind  is  oriented  in  the 
positive  earth  axis  direction. 
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Total  wind  velocity  is  composed  of  the  ambient  wind  and  the  non-random 
disturbances  introduced  by  the  airwake.  This  section  defines  the  non-random 
wind  components  in  aircraft  body  axes. 

VOLAR  is  currently  set-up  to  initialize  on  airspeeds  instead  of  inertial 
speeds.  Equations  for  initializing  total  inertial  velocity  components  are 
contained  in  this  section.  They  are  calculated  once  only.  This  calculation 
is  controlled  by  the  ICEES  parameter. 

"Airspeed" 

Wind  velocity  components  are  subtracted  from  inertial  velocity  components 
to  yield  airspeed  components  in  the  aircraft  body  axis  system.  Angle-of-attack, 
sideslip,  and  total  airspeed  are  calculated  according  to  standard  conventions. 

"Look  Up  Aero,  and  Open  Loop  Parameters" 

Aerodynamic  stability  and  control  derivatives  are  stored  in  tables  as  a 
function  of  airspeed.  This  section  of  the  program  performs  a table- look-up 
of  these  derivatives  at  the  current  value  of  airspeed.  Subroutine  TLU 
evaluates  several  dependent  variables  with  one  call.  Fifteen  parameters,  YV 
through  NDR  (R  array,  segment  10),  are  determined  by  the  single  call  to  TLU. 
Open  loop  aileron  and  rudder  control  deflections  are  set  to  zero. 

"Range  From  Ship" 

Range  from  ship  is  arbitrarily  defined  here  as  the  distance  between  the 
x earth  axis  positions  of  the  ship  and  airplane.  "Range",  as  defined  here, 
provides  some  measure  of  hov-  close  the  airplane  is  to  the  ship. 

"Lateral  Position  Error" 

The  aircraft's  commended  lateral  position  is  the  lateral  position  of 
the  ship.  Error  in  position  is  the  difference  between  the  two. 

This  section  begins  the  build-up  of  parameters  associated  witn  the 
calculation  of  pilot  control  inputs.  Figure  8 (Volume  1,  page  40)  is  a 
schematic  of  the  classical  pilot  model  employed  for  the  AV-8A  at  low  speed. 

"Roll  Attitude  Error" 

Corrmanded  bank  angle  is  a gain  multiplied  by  lateral  position  error. 

Bank  angle  error  Is  the  difference  between  the  corrmanded  and  actual  hank 
angles.  This  error  is  modified  by  a bias  term  to  account  for  the  initial 
trim  bank  angle.  A non-zero  initial  bank  angle  is  necessary  to  trim  out  the 
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side  forces  on  the  airplane.  These  side  forces  result  from  the  steady  side 
velocity  required  for  tracking  the  ship. 

"Heading  Error" 

The  commanded  aircraft  heading  is  zero  degrees.  Heeding  error  is  the 
difference  between  commanded  and  actual  heading. 

"Roll  Control" 

The  pilot's  roll  control  input  is  calculated  here.  From  Figure  8 
(Volume  I,  page  40)  the  appropriate  transfer  function  is: 

Yp  _ -112.0  (s  + 0.1)  (s  - 10.0) 

“ (s  +10.0)  (s  + 10.0) 

This  transfer  function  was  obtained  from  the  following  transfer  function. 

v (T|i  ^ + 1-0)  (2/t  . - s) 

% " K1  (Ta  s + 1.0)  f2/Td  + s) 

The  quantity  Kj  (tl^  s + 1)  is  the  compensation  applied  by  the  pilot.  A lag, 
(t  s + 1),  is  introduced  to  represent  actuator  and  neuromuscular  lags.  The 

^ I 

final  numerator  and  denominator  factors  provide  a first  order  Pade  approxima- 
tion to  a pure  time  delay.  This  delay  represents  the  time  it  takes  the  pilot 
to  generate  an  output  strategy  following  sensing  of  the  error.  Time  delays 
are  a standard  feature  in  most  pilot  models. 

As  presented  in  block  diagram  form  below,  the  pilots'  output,  DAP,  must 
pass  through  a control  system  limiter. 


Subroutine  DESCRIB  calculates  the  non-linear  describing  function  which 
represents  the  limiter.  A necessary  input  to  DESCRIB  is  the  standard  devia- 
tion of  the  signal  DAP.  A general  equation  for  calculating  the  standard 
deviation  of  a parameter  which  is  a function  of  the  state  variables  but  not 
a state  variable  itself  is: 


y 


where  y = f(xj,  x2,  ....  xR) 

and  P is  the  Covariance  matrix  associated  with  the  state  vector  x^ 

Normally  the  standard  deviation  of  the  control  would  be  evaluated  by 
describing  functions  and  the  mean  would  be  determined  by  applying  the  true 
limiter.  This  was  not  done  here.  For  this  application,  the  difference  in 
results  was  insignificant  and  so  the  modified  mean  (i.e.  its  describing 
function  representation)  was  used  to  reduce  the  number  of  calculations. 

"Yaw  Control" 

This  section  is  so  similar  in  layout  to  the  previous  one  it  will  not  be 
discussed.  The  output  is,  of  course,  the  yaw  control  from  the  pilot. 


"Ship" 

Only  the  mean  lateral  ship  motion  is  considered  for  this  example.  The 
shiD  is  steaming  with  velocity  V5  in  a direction  PS  1 5 with  respect  to  north. 
When  integrated,  the  differential  equation  in  this  section  yeilds  the  y earth- 
axis  position  of  the  ship.  The  x and  z earth-axis  positions  were  defined  in 
previous  sections  of  this  subroutine;  one  by  table  look-up  end  the  other  by 
equation. 


"Airframe" 

This  section  of  coding  defines  the  differential  equations  which  describe 
the  aircraft's  mean  trajectory.  Complete  non-linear  equations  are  used.  This 
creates  a small  problem  because  the  aerodynamic  data  are  furnished  as  dimensional 
stability  derivatives  which  give  perturbation  forces  and  moments.  Total  forces 
and  moments  are  required  to  drive  the  mean  equations  of  motion.  To  demonstrate 
how  this  is  overcome,  consider  the  V equation. 

Some  definitions: 

Let  Y_  be  the  total  side  force  divided  by  mass  at  some  reference  time, 
o 


and  SR  are  the  pilot  commanded  positions 

6.  and  <5d  are  open  loop  inputs 
m0L  kql 


v„  Is  the  perturbation  airspeed  along  the  y body  axis 

aS 
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p and  r are  the  perturbation  roll  and  yaw  rates. 
P and  R are  the  total  roll  and  yaw  rates 
U,  V,  and  W are  total  airspeed  components 


Based  on  these  definitions  proceed  as  follows: 

1.  the  total  control  deflection  is  given  as: 
S4  ■ «A  * V 
<5r  = <So  + 6n 


the  incremental  or  perturbation  force-over-mass  is: 


AY  = 


Y v +Yp  + Yr 
v as  pK  r 


3.  the  total  force-over-mass  on  the  airplane  is: 

Y = Yq  + AY  + Y6a  64  + Y6r  65 

4.  the  rate  of  change  of  side  velocity  is: 

V = Y + g cos(0)  si n (<J>)  - R U + P W 

Thus  the  total  y-axis  velocity  component,  V,  can  be  found  by  integrating  V. 
"Update  Sub  Zero  Values" 

The  zero  subscripted  values  define  the  reference  trajectory  about  which 
quasi  linearisation  occurs.  Perturbation  quantities  are  found  by  subtracting 
the  reference  value  from  the  total  value.  The  reference  values  are  updated 
prior  to  each  integration  increment  and  remain  fixed  during  the  increment. 

Zero  subscripted  terms  also  appear  as  elements,  or  parts  of  elements,  in  the 
F matrix.  (F  is  the  system  matrix  appearing  in  the  covariance  propagation 
equation. ) 

PROP 

Subroutine  PROP  has  only  one  purpose;  to  furnish  the  integration  sub- 
routine with  P.  The  prospective  user  should  keep  this  in  mind  when  reading 
through  the  subroutine.  Familiarity  with  state  variable  notation  Is  helpful. 

As  an  overview,  consider  the  system  governed  by  the  differential  equation, 
z * F z_  + w 

where  z is  the  state  vector, 

F is  the  system  matrix 

GM  Is  a control  or  coupling  matrix 

w is  white  noise 
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The  covariance  matrix,  P,  is  defined  by: 

P = E { <2  - Z)  (2  - I)' ) 

and  obtained  by  intenrating  the  following  differential  equation. 

P = F P + PFT  + GM  Qm  gJ 
A T 

Qm  = E{ww> 

Each  of  the  titles  preceding  a group  of  calculations  in  subroutine  PROP  is 
explained  below,  with  the  titles  repeated  in  quotation  marks. 

"Miscellaneous  Calculations" 

This  group  of  calculations  provides  some  frequently  used  parameters. 

They  are  grouped  here  for  "one  time  per  call"  calculation.  This  minimizes 
calls  to  the  trigonometric  functions,  etc. 

"FI  Matrix" 

The  FI  matrix  is  the  aircraft  system  matrix.  Aircraft  states  are  v,  p, 
r,  <J>,  \Ji,  and  y.  FI  is  developed  from  the  standard  linearized  equations  of 
motion.  Refer  to  Section  2.0  in  Volume  I for  a description  of  the  linearized 
aircraft  equations. 

"G1  Matrix" 

G1  is  the  control  matrix  for  the  airplane.  The  airplane  is  represented 
by: 

x = Fj  x + Gj  u^ 

where  u is  a vector  of  control  inputs. 

"Gamma  Matrix" 

Gamma  is  the  matrix  for  coupling  the  colored  noise  disturbances  with  the 
airplane.  If  w'  represents  a vector  of  colored  noises,  the  equations  describing 
the  airplane  motion  are. 

x = Fj  x ♦ Gj  u + r w' 

"A  Matrix" 

The  A matrix  and  the  B matrix  (next  section)  define  the  noise  coloring 
processes.  In  this  example,  noise  is  composed  of  Dryden  atmospheric  turbulence, 
plus  the  ship's  airwake.  Equations  for  the  noise  take  the  form: 


— 


U 

0 

11 

0 

0 

D 


- 


fl 


D 


-i 
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w'  = Aw'  + B w 

where  w'  is  the  colored  noise  (output)  and  w is  white  noise  (input).  Dryden 
turbulence  is  represented  by  a first  order  filter  (in  each  axis).  Each  axis 
of  ship  airwake  requires  a second  order  filter. 

11 B Matrix" 

The  B matrix  is  explaned  in  the  preceeding  section. 

"C  Matrix" 

C and  Rl  are  used  to  define  the  control  input,  £. 
jj  = Cx.  + Rl  u. 


x and  u^  are  explained  above. 


11 F Matrix11 

F is  the  system  matrix  for  the  augmented  vector 


F = 


0 

C 


r 

A 

o 


( x. , w ' , u) : 


The  subroutine  LAYIN  is  designed  to  store  submatrices  (Fj,  F,  etc.)  into  a 
larger  matrix  (F).  3y  building  F from  its  submatrices;  (1)  the  user  has 
more  visibility  of  the  pieces,  (2)  there  is  less  chance  of  an  error  in  row 
or  column  placement  of  individual  elements,  (3)  elemental  model  additions  or 
size  changes  can  be  Incorporated  quickly. 

"GH  Matrix” 

Gm  Is  the  multiplier  of  white  noise  Inputs.  For  this  example  G(i,  = B. 

"QM  Matrix” 

Qm  Is  defined  by: 

Qm  = E (w  wT  } 

The  (1,  1)  element  corresponds  to  the  Dryden  turbulence  end  the  remaining 
diagonal  elements  are  associated  with  the  ship's  airwake. 
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“Propagate  Covariances" 

These  operations  develop  P from  F,  GM,  and  Q^.  Matrix  operations  are 
identified  in  subroutine  MATRIX.  MATRIX  is  written  in  COMPASS;  but  FORTRAN 
equivalent  subroutines  appear  as  comments  proceeding  each  ENTRY  point. 


APPENDIX  B 


LISTING  FOR  CLASSICAL  PILOT  MODEL 
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D.l  Decomposition  Into  Odd  and  Even  Components 

The  general  nonlinearity  is  represented  by  N straight-line  segments. 
The  number  of  segments  in  the  odd  component  is  also  N,  and  there  are  also  N 
segments  in  the  even  component. 

Denoting  the  general  nonlinearity  by  v = f(x):  then  the  odd  and  even 
components  are  respectively  aiven  b.v 


yQ(x)  = 1/2  y(x)  - v(-x) 
yE(x)  = 1/2  y (x)  + v(-x) 


Note  that  for  x = 0 this  gives  yo(0)  = 0:  vE(0)  = v(0).  It  is  convenient  to 
separate  out  the  offset  C,  as  indicated  bv  the  following  equations 

V = yo(x)  = fo'(x) 

yE'(x)  = yE(x)  - C = fE'(x) 
where  C = y(0),  as  shown  in  Figure  D.l. 
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Figure  0.2  Separation  of  the  Nonlinearity  Into  Fven,  Odd,  and  Bias  Components 
D. 2 Approximation  to  the  Odd  Component 


The  input  is  assumed  to  consist  of  a bias,  B,  plus  a random  Gaussian 
component  of  variance  a2,  denoted  by  xR.  The  aeneral  form  of  the  odd  component 
— is  shown  in  Figure  A. 5. 3. 

x(t)  = B + xR(t) 

The  best  least-sauares  approximation  to  yQ'  is 
yoA  = NB  B + NR  xR^t^ 

where  NR  and  NR  are  given  by  the  expressions  overleaf  (from  p.  579  of 
Reference  D.l). 

0 
0 
L 
LI 

0 
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6,  + B 


6,  - B 


Nr  - 2e  PE(|)  * 2mp  - mj  * ^ (mf  - i»h1)[p1  (-^-J  * PII-^-)] 


6.  + B 


[2PI(|)-l]  * 2mp  - mj  * |T  («,  - m1+1)[s<^)  - P (^)] 


B ” B 


where  mi  is  the  slope  of  the  i-  segment  (See  Figure  D.3) 


PF,  PI;  and  G are  defined  below 


PF(w)  = ( — T ] ■ 

Ztt 


-l/2w2 


PI (w)  = /w  PF  (n)  dn 


G(w)  = wPI(w)  + PF(w) 

Note  that  D = 0,  because  we  have  separated  out  the  bias  C. 

D.3  Approximation  to  the  Even  Component 

The  best  least-souares  approximation  to  y^'  is  obtained  bv  subtracting 
two  describing  functions  given  on  p 585  and  586  of  Reference  0.1.  The  results 
of  subtracting  describing  function  No.  55  from  describing  function  No.  52 
reduce  to: 


6,  + B 


6,  - B 


Nn(j)  = m.  PI  (-J— -)  - m.  PI  ) 


5.  6.+B  m.  a 6.  - B 

nbU)  ■ -Zmj  i * »,  | P (J^_) 

where  mj,  6j  are  defined  as  illustrated  in  Figure  D.4,  for  a non-qeneral 
even  characteristic  with  only  one  breakpoint  on  each  side  of  the  v^  axis. 

D.l.  Gelb,  A.,  W.  E.  Vander  Velde,  Multiple-Input  Descrioing  Functions  and 
Nonlinear  System  Design,  McGraw-Hill,  1968. 


:•  j 


AD-A066  173  VOUGHT  CORP  DALLAS  TEX  F/G  1/2 

volar:  a digital  computer  program  for  simulating  vstol  aircraft— ETC <U) 

DEC  78  J WOLKOVITCH*  R B BRASSELL  N62269-77-R-0389 

UNCLASSIFIED  2-32000/8R-51672-VOL-2  NADC-77123-30-VOL-2  NL 


4 *4 

ADA 

066)73 

El 

H 

1 

M 

HI 

H 

p 

p 

n 

Hi 

9 

1 

1 

fj 

pi 

H 

H 

P 

H 

9 

H 

ij 

gg 

If 

■ 

H 

f| 

9 

n 

END 

DATE 

FILMED 

5 79 

DOC 

I 


Figure  D.4  Non-General  Even  Characteristic 

In  Figure  D.4  is  the  only  breakpoint.  For  the  general  even  component 
there  will  be  0 breakpoints  as  shown  in  Figure  A. 5,  hence  the  overall  apDroxi 
mation  is: 


where: 

j=Q 

V - r NR(j) 

R J=1  R 

and 

j=Q 

V = E No(j) 

B j=l  B 

Hence  the  total  output  is  approximated  by: 

- B <N,'  * Nb)  * *R(t)(NR-  * N„)  ♦ C 

D.4  Routines  for  Computing  Describing  Functions  of  General  Single-Valued 
Nonlinearities  Represented  by  Straight-Line  Segments 

Three  FORTRAN  routines,  TABRD,  DESCRIB  and  ERF  provide  a capability  for 

generating  describing  functions  for  general  single-valued  nonlinearities 

which  can  be  expressed  as  functions  of  as  many  as  three  independent  variables.* 

One  COMMON  block  is  required  to  transfer  data  between  routines.  The  procedures 

for  invoking  the  proper  routines  and  initializing  the  COMMON  block  are 

described  below. 

Method 

For  ease  of  reference  the  key  formulas  of  the  preceding  section  are 
repeated  below. 

Denoting  the  general  nonlinearity  by 

V = F(X) 

then  the  odd  and  even  components  are  given  bv 

Y0(X)  = 1/2 (Y ( X)  - Y(X) ) 

Ye(X)  = 1/2 ( Y (X ) + Y ( -X ) ) 

*For  VOLAR  purposes  only  one  independent  variable  is  assumed. 
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for  X = 0,  this  yields 


Yo(0)  = 0 
Ye(0)  = Y(0) 


It  is  convenient  to  separate  out  the  offset  yielding 


0 

l n 


Y*(X)  = Yo(X) 

Yq(X)  = Ye(X)  - Y(0) 

Assuming  the  input  consists  of  a bias,  B,  and  a random  Gaussian  component 
of  variance  o2,  denoted  by  Xr>  the  best  linear  least  squares  approximation  to 

*o  is: 

y0  ■ Nb  B * Nr  Xr<‘> 


where: 


6,  + B 


6.  + B 


6,  - B 


" Mi  * I <Hi  - Hm>  [g<-4->  - - g(-4-i] 


6.  + B 


6,  - P 


r-i  r o.  t D r 1 

- Mi + { <Mi  - Vi>  [PI<-4-> ♦ PI<-4->] 


Mi  = Slope  of  the  ith  odd  component  segment 


P = Index  of  last  odd  component  linear  segment 


G(VI)  = WPI(W)  + PF(W) 
PI (W)  = "/  PF(n)  dn 


PF<“>  * V2 ~ e 


Max  value  of  the  first  independent  variable 
associated  with  the  ith  odd  component  segment 


NOTE:  Only  segments  where  6^0  are  considered. 

The  best  linear  least  souares  approximation  to  is: 


N£  B + n;  xr(t) 


where: 


6.  + B 


6.  - B 


^ - «,>  [pi<-V->  - pi<-4->] 

P-1  r-26.  6.  + B 6,  - B 

Ul  ("hi  - ", ) [nr  + ff(R(-V-)  * s(-V-»] 


Mj  - Slope  of  the  ith  even  component  segment 

P = Index  of  last  even  component  linear  segment 

G,  PI,  PF  as  defined  above 

6..  = Min  value  of  the  first  independent  variable 

associated  with  the  ith  even  component  segment 

NOTE:  Only  segments  where  6^>0  are  considered. 

The  approximated  output  bias  is  given  by 


Ya  s B(Nb  + Nb>  + Y(0) 

The  approximated  output  variance  is  given  by 

°2  8 (Nr  ♦ n;)2 


D.5  Common  Block  Initialization 

It  can  be  seen  from  the  discussion  on  methodology  above  that  the  data 
that  must  be  preserved  includes: 

(1)  Values  of  the  second  and  third  independent 
variables  for  interpolation  purposes.* 

(2)  A subset  of  the  first  indepent  variable. 

When  generating  the  even  and  odd  components 
of  the  nonlinearity,  the  x values  are 
symetrically  disposed  around  x = 0.  Once 
generated,  only  the  disposed  values  of 

x i 0 need  be  kept. 

(3)  Slopes  of  the  even  and  odd  components  corresponding 
to  the  stored  values  of  x. 

(4)  The  offset  that  was  separated  out  of  the  even 
components . 

The  method  of  data  storage  within  the  COMMON  block  is  discussed  below. 

C0MM0N/GDF/NTABLE , INDEX(m) .TABLE (n ) 

NTABLE  - Number  of  tables  that  were  read. 

INDEX  - Array  holding  pointers  into  TABLE  (must  be 
dimensioned  at  least  3*NTABLE+1). 

INDEX(3*(I-1)+1)  - Index  into  TABLE  where  data  for 
table  I begins. 

INDEX(3*(I-l)+2)  - Number  of  values  of  2nd  independent 
variables  in  table  I. 

INDEX( 3*1 ) - Number  of  saved  values  of  x In  table  I. 


*For  VOLAR  purposes  only  one  independent  variable  is  used.  Zeros  were 
input  in  the  table  for  the  second  and  third  independent  variables. 


x3  VALUE  (1  ENTRY) 


y 

E 

I 

0 

0 

I 

I 

II 
D 

E 

I 

I 


X2  VALUES 

(INDEX(3*(I-l)+2)  ENTRIES) 


SAVED  Xj VALUES 
( INDEX (3* I ) ENTRIES) 


OFFSET  SEPARATED  FROM  EVEN 
COMPONENTS  ( INDEX(3*(I-l)+2)  ENTRIES) 


ODD  COMPONENT  SLOPES  FOR 
1st  Y VALUE ( INDEX (3*1 )-l  ENTRIES) 


EVEN  COMPONENT  SLOPES  FOR 
1st  Y VALUE ( INDEX( 3*1  )-l  ENTRIES) 


c 

o 

ODD  SLOPES  FOR  INDEX(3*(l-l)+2)th 
Y VALUE  ( I NDE X( 3*1 ) -1  ENTRIES) 

EVEN  SLOPES  FOR  INDEX(3*(I-i)+2)th 
Y VALUE ( I NDEX ( 3*1 )-l  ENTRIES) 


TABLE  must  be  dimensioned  large  enough  to  hold  required  data  generated 
from  all  input  tables. 
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D.6  Routine  Description 
TABRD 

This  routine  is  responsible  for  calculating  the  even  and  odd  components 
of  the  nonlinearity,  calculating  the  component  slopes  and  storing  the  data 
into  the  COMMON  block.  The  calling  statement  is: 

CALL  TABRD(X,Y,NPTS) 

where: 

X - Array  of  independent  variables. 

Y - Array  of  dependent  variables. 

NPTS  - Number  of  elements  in  the  X and  Y arrays. 

DESCRIB 

This  routine  is  responsible  for  calculating  elements  of  the  describing 
function.  The  calling  statement  is: 

CALL  DESCRIB( ITAB, BIAS, SIGMA, 0BIAS,0VAR,REPR0) 

where: 

ITAB  = Table  number  for  which  describing  function 
is  to  be  generated. 

BIAS  = Input  BIAS. 

SIGMA  = Input  SIGMA. 

0BIAS  = Output  BIAS. 

0VAR  = Output  variance. 

REPR0  = Gain  on  output  sigma. 
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The  Input,  ALPHA,  consisted  of  a bias  of  10  deg.  and  a variance  of  1 deg. 
The  describing  function  yielded  a mean  for  the  output,  CM,  of  -.139988  and  a 
variance  of  0.000052097.  An  alternate  method  of  Inputting  1000  Gaussian- 
distributed  random  numbers  with  a mean  of  10  and  a variance  of  1 was  checked 


yielding  an  output  mean  for  CM  of  -.14  and  a variance  of  0.00005279. 


D.l.  Gelb,  A.,  anu  A.  E.  Varuler  Velde,  Multi ple-Inout  Describing  Functions 
and  Nonlinear  System  Design,  McGraw-Hill,  1968. 

D.2.  Nave,  R.  L.,  "Progress  Toward  a Computerized  VSTOL/Small  Platform 
Landing  Dynamics  Investigation  Model,"  NADC  Report  77D24-30, 

1977. 


E.l  The  I Array 


I ARRAY  VARIABLE 

SEGMENT  ELEMENT  NAME 


DESCRIPTION 


1 


1 

1 

1 

1 

1 

1 

1 

1 


1 


1 

ICEES 

An  internal  'flag'  to  keep  track  of  what 
stage  the  integration  is  in  and  to  identify 
the  very  first  pass  through  the  program 
(ICEES  = 1). 

n 

C. 

NEQU 

Total  number  of  differential  equations  to 
integrate. 

3 

NSTATES 

Number  of  mean  equations  to  integrate. 

4 

MPA 

Order  of  the  P matrix  (NPA  by  NPA). 

5 

NQM 

Order  of  the  QM  matrix  (NO!'1  by  NQM). 

6 

NGM 

Row  dimension  of  GM  matrix  (NGM  by  NOM). 

7 

NPE 

Order  of  the  PE  matrix  (NPE  by  NPE). 

8 

NPX 

Order  of  tne  PXHAT  matrix  (NPX  by  NPX). 

9 

N0AW 

User  specified  control  of  airwake  calcula- 
tions . 

N0AW  - 0,  calculate  airwake 

N0AV!  = 1,  no  airwake  calculations  desired 

10 

N0 : N IT 

User  specified  control  of  program  flow 

N0INIT  = -1,  fix  means  and  integrate 
covariance  matrices  only 

N0INIT  - 0,  call  SOLVR  to  calculate 

steady  state  value  of  P at  the 
conditions  defined  by  means 

N0INIT  = i,  no  special  initialization 
calculations  desired  - run  a 
normal  time  history 
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II 

I E 

11 

G 

hi 


Li 

r 

| E 

h 


I 

ARRAY 

VARIABLE 

DESCRIPTION 

SEGMENT 

ELEMENT 

NAME 

4 

1-14 

INDX1 

R array  locations  of  elements  to  be 
integrated.  Setup  by  Option  6. 

5 

1-14 

INDX2 

R array  locations  of  elements  that  are 

outputs  of  integrators.  Setup  by  Option  6. 


T 

, 

SEGMENT 

ARRAY 

ELEMENT 

VARIABLE 

NAME 

TT 

ii 

1 

DESCRIPTION 

6 

1 

I.NDXPO 

R array  location  of  the  first  element  of 

the  PD0T  matrix.  Setup  by  Option  t.  ■* 

6 

2 

INDXP 

R array  location  of  the  first  element  of 

the  P matrix.  Setup  by  Option  6.  11 

6 

3 

INDXPED 

R array  location  of  the  first  element  of 
the  PED  matrix.  Setup  by  Option  6. 

6 

4 

INDXPE 

R array  location  of  the  first  element  of 
the  PE  matrix.  Setup  by  Option  6. 

6 

5 

INDXP XD 

R array  lecatior  of  the  first  element  of 
the  PXHATD  .matrix.  Setup  by  Option  6. 

6 

6 

INDXPX 

1 . 

R array  location  of  tie  first  element  of 
the  PXhAT  .liatrix.  Setup  by  Option  6. 

U 

D 

D 

D 

W 


L 

0 

0 

0 

[] 

y 

E 

0 

P 

f 

D 

D 
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E.2  The  R Array 


R 

SEGMENT 

ARRAY 

ELEMENT 

VARIABLE 

NAME 

UNITS 

DESCRIPTION 

l 

1 

TIME 

t 

T ime 

1 

2 

TO 

t 

Initial  time;  time  for  which 
reference  values  are  defined 

1 

3 

DTI  ME 

t 

Integration  interval 

1 

4 

DT0UT 

t 

Output  file  written  every 

D70UT  seconds 

1 

5 

DT0UT2 

t 

Print  interval  for  MISCAL 
output 

1 

6 

TMAX 

t 

Upper  integration  limit 

i 
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R 

SEGMENT 

ARRAY 

ELEMENT 

VARIA6LE 

NAME 

UNITS 

DESCRIPTION 

3 

1 

UD 

1/t 

Inertial  speed  along  airplane 

X body  axis 

3 

2 

WB 

1/t 

Inertial  speed  along  airplane 

Z tody  axis 

3 

3 

QB 

rad/t 

Airplane  pitch  rate 

3 

4 

THETA 

rad 

Airplane  pitch  attituae 

3 

t 

XAPP 

1 

Airplane  position,  X earth  axis 

3 

6 

ZAPP 

1 

Airplane  position,  Z earth  axis 

3 

7 

DFRPM 

rpm 

Engine  rpiF 

3 

8 

XSP 

1 

Strip  position,  X earth  axis 

3 

9 

THETAS 

rad 

Ship  pitch  attitude 

3 

10 

VB 

1/t 

Inertial  speed  alone  airplane 

Y body  axis 

3 

11 

PB 

rad/t 

Airplane  roll  rate 

3 

12 

RB 

rad/t 

Ai rplane  yaw  rate 

3 

13 

YAPP 

1 

Airplane  position,  Y earth  axis 

3 

14 

PHI 

rad 

Airplane  roll  attitude 

3 

15 

PSI 

rad 

Airplane  heading 

3 

16 

YSP 

1 

Ship  position,  Y earth  axis 

R ARRAY 

SEGMENT  ELEMENT 


VARIABLE 

NAME 


UNITS 


DESCRIPTION 


l 


D 

1 

UD0T 

1/t2 

Airplane  acceleration  along  X 
body  axis 

2 

u 

2 

WD0T 

1/t2 

Airplane  acceleration  along  Z 
body  axis 

2 

| 

3 

QD0T 

rad/t2 

Airplane  pitch  acceleration 

2 

4 

TD0T 

rad/t 

Airplane  pitch  attitude 
rate-of-change 

L 2 

5 

XAPPD0T 

1/t 

Airplane  velocity  along  X 
earth  axis 

h 

6 

ZAPPD0T 

1/t 

Airplane  velocity  along  Z 
earth  axis 

fi 

7 

DFRPMD 

rpm/t 

Engine  rpm  rate-of-change 

2 

8 

XSPD0T 

1/t 

Ship  velocity  along  X earth  axis 

D 2 

9 

ZSPD0T 

1/t 

Ship  velocity  along  Z earth  axis 

r 

10 

TSPD0T 

rad/ 1 

Ship  pitch  attitude  rate-of-change 

R ARRAY 

SEGMENT  ELEMENT 

VARIABLE 

NAME 

UNITS 

DESCRIPTION 

4 

1 

UO 

1/t 

Reference  value  of  UB  used  for 
local  linearization 

4 

2 

WO 

1/t 

Reference  value  of  >/B  used  for 
local  linearization 

4 

3 

QO 

rad/t 

Reference  value  of  QB  used  for 
local  linearization 

4 

4 

THETAO 

rad 

Reference  value  of  THETA  used 
for  local  linearization 

4 

5 

XAPPO 

1 

Reference  value  of  XAPP  used  for 
local  linearization 

4 

6 

ZAPPO 

I 

Reference  value  of  ZAPP  used  for 
local  linearization 

4 

7 

XO 

1/t2 

lFx/rnass 

4 

8 

ZO 

1/t2 

rFz/mass 

4 

9 

MO 

rad/t2 

EMy/ I y 

4 

10 

UBASO 

1/t 

Reference  value  of  DBAS 

4 

11 

WBASO 

1/t 

Reference  value  of  '/BAS 

4 

12 

WASDO 

1/t2 

Reference  value  of  WASD 

4 

13 

THETAAO 

rad 

Reference  value  of  THETA, 

4 

14 

PS  I AO 

rad 

Reference  value  of  PSIA 

4 

15 

PHIAO 

rad 

Reference  value  of  PHIA 

4 

16 

RCDO 

1/t 

Reference  value  of  RCD 

4 

17 

VO 

1/t 

Reference  value  of  VB  used  fcr 
local  1 inearization 

4 

18 

PO 

rad/t 

Reference  value  cf  PB  used  for 
local  linearization 

4 

IS 

RO 

rad/t 

Reference  value  of  RB  used  for 
local  linearization 

E-8 
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R ARRAY 
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VARIABLE 

NAME 


UNITS 


GESCRIPTIGN 


4 

20 

YAPPO 

1 

Reference 

value  of  YAPP 

4 

21 

PHIO 

rad 

Reference 

value  of  PHI 

4 

22 

PSIO 

rad 

Reference 

value  of  PSI 

4 

23 

YO 

1/t2 

LFy/mass 

4 

24 

LO 

rad/t2 

4 

25 

NO 

rad/t2 

EM‘/IZ 

4 

26 

VBASO 

1/t 

Reference 

vo.lue  of  VBAS 

L 

li 

I 
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R ARRAY 

SEGMENT  ELEMENT 

VARIABLE 

NAME 

UNITS 

5 

1 

SUMFXO 

f 

5 

2 

SUMFZO 

f 

5 

3 

SUMMO 

fl 

5 

4 

SUMFX 

1/t2 

5 

5 

DFX 

1/t2 

5 

6 

SUMFZ 

1/t2 

5 

7 

DFZ 

1/t2 

5 

8 

SUMM 

rad/t2 

5 

9 

0M 

rad/t2 

5 

10 

SUMFYO 

f 

5 

11 

SUMLO 

fl 

E- 10 


DESCRIPTION 


l 

I 

1 


Initial  summation  of  forces 
along  the  airplane  body  axis. 
Required  input  if  means  are 
propagated. 

Initial  summation  of  forces 
along  the  airplane  body  axis. 
Required  input  if  means  are 
propagated. 

Initial  summation  of  moments 
about  the  airplane  body  axis. 
Required  input  if  means  are 
propagated. 

Summation  of  'forces'  along 
airplane  X body  axis. 


Incremental  'force'  due  to 
airplane  departure  from 
reference  trajectory. 

Summation  of  'forces'  along 
airplane  Z body  axis. 

Incremental  'force'  due  to 
airplane  departure  from 
reference  trajectory. 

Summation  of  'moments'  about 
airplane  Y body  axis. 

Incremental  pitching  'moment* 
due  to  airplane  departure  from 
reference  trajectory. 

Initial  summation  of  forces 
along  the  airplane  Y body  axis. 
Required  input  if  means  are 
propagated. 

Initial  summation  of  moments 
about  the  airplane  X body  axis. 
Required  input  if  means  are 
propagated. 


I 

B 

0 


1 


L 

U 

• m 

0 

D 

11 
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SEGMENT 

ELEMENT 

NAME 

UNITS 

DESCRIPTION 

5 

12 

SUMNO 

fl 

Initial  summation  of  moments 
about  the  airplane  Z body  axis 
Required  input  if  means  are 
propagated. 

5 

13 

SUMFY 

1/t2 

Summation  of  'forces'  along 
airplane  Y body  axis. 

5 

14 

DFY 

1/t2 

Incremental  'force'  due  to 
airplane  departure  from 
reference  trajectory. 

5 

15 

SUML 

rad/t2 

Summation  of  'moments'  that 
produce  angular  acceleration 
only  about  the  airplane  X body 
axis. 

5 

16 

DL 

rad/t2 

Incremental  'moment'  due  to 
airplane  departure  from 
reference  trajectory. 

5 

17 

SUMN 

rad/t2 

Summation  of  'moments'  that 
produce  angular  acceleration 
only  about  the  airplane  Z body 
axis . 

5 

18 

DN 

rad/t2 

Incremental  'moment'  due  to 
airplane  departure  from 
reference  trajectory. 

1 


li 

0 

I 


E-ll 


I 


R ARRAY 

SEGMENT  ELEMENT 

VARIABLE 

NAME 

UNITS 

DESCRIPTION 

6 

1 

WEIGHT 

ml/t2 

I 

Airplane  weight 

6 

2 

G 

1/t2 

Gravitational  constant 

1 

6 

3 

IX 

ml  2 

Airplane  body  axis  inertia 
component 

6 

4 

IY 

ml 2 

Airplane  body  axis  inertia 
component 

6 

5 

IZ 

ml  2 

. 

Airplane  body  axis  inertia 
component 

6 

6 

IXZ 

ml 2 

Product  of  inertia 

6 

7 

LM 

1 

Overall  length  of  ship 

6 

8 

WM 

1 

Maximum  beam  of  ship 

6 

9 

LTDM 

1 

Distance  from  bow  to  landing 
pad  center  measured  along  ship 
centerline 

6 

10 

LS 

1 

Scaled  ship  length 

6 

11 

LTDS 

1 

A scaled  LTDM 

I 1 

I (1 

D 

0 
0 

0 

0 
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R ARRAY 
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VARIABLE 

NAME 


UNITS 


DESCRIPTION 


1 

0MEGAU 

WHITE 

NOISE 

2 

0MEGAW 

3 

KU 

WHITE 

NOISE 

4 

KW 

5 

KFRPM 

6 

KCDN 

7 

TAUENG 

t 

8 

TAUCDN 

t 

9 

TAUD 

t 

10 

KDEU 

1/1  ' 

11 

KDEW 

1/1 

12 

KDEQ 

1/rad 

13 

KDET 

1 

t rad 

14 

KDEX 

1 

FT 

15 

KDEZ 

1 

t 1 

16 

KDEN 

1 

t rpm 

17 

KDEDE 

1 

t 

18 

KDEDT 

1 

t 

KU 

S + 0MEGAU 


KW 

S + 0MEGAW 


DRYDEN  TURBU- 
LENCE ALONG  AIR- 
PLANE X BODY  AXIS 


DRYDEN  TURBU- 
LENCE ALONG  AIR- 
PLANE Z BODY  AXIS 


Gain  on  input  to  engine  model 
Gain  on  input  to  nozzle  actuator 
Engine  time  constant 
Nozzle  actuator  time  constant 
Nozzle  actuator  time  constant 


Longitudinal  optimal  pilot 
feedback  gains  on  states 
into  pitch  control 


[ 

D 

Q 


31 

32 

33 

34 

35 

36 

37 


K1 

K2 

KY 

B1 

B2 

TL1 

TL2 


Cont'd 

R 

SEGMENT 

ARRAY 

ELEMENT 

VARIABLE 

NAME 

UNITS 

7 

19 

KDTU 

N 

1/1 

7 

20 

KDTW 

1/1 

7 

21 

KDTQ 

1/rad 

7 

22 

KDTT 

1 

t rad 

7 

23 

KDTX 

1 

t 1 

7 

24 

KDTZ 

1 

t 1 

7 

25 

KDTN 

1 

t rpm 

7 

26 

KDTDE 

1 

t 

7 

27 

KDTDT 

1/t 

7 

28 

0MEGAV 

WHITE 

NOISE 

7 

29 

30 

KV 

TAUA 

7 

t 

DESCRIPTION 


Longitudinal  optimal  pilot 
feedback  gains  on  states 
into  thrust  control 


K V 

S + 0MEGAV 


DRYDEN  TURBU- 
LENCE ALONG  AIR- 
PLANE Y BODY  AXIS 


Actuator  time  constant 
ROLL  CONTROL: 


ERPHI 


K1 

rTLl*S+Bl 

( -S+2/TAUD) 

. 

[TAU*S+1) 

S+2/TAUD) 

DAP 


YAW  CONTROL: 
ERPSI 


K2(TL2*S+B2)(-S+2/TAUDl 
i TAU*S+1 ) ( S+2/TAUD) 


DRP 


1 

1 

I 

I 

I 


i: 


i . 


I 

l 

1 

fS 


o 

D 

E 


E- 14 


71 

11 

R 

ARRAY 

VARIABLE 

UNITS 

DESCRIPTION 

B 

SEGMENT 

ELEMENT 

NAME 

I n 

8 

1 

ALT 

1 

Altitude 

! " 

8 

2 

VAIR 

1/t 

Ai rspeed 

D 

8 

3 

UBAS 

1/t 

Airplane  X body  axis  component 
of  airspeed 

8 

4 

VBAS 

1/t 

Airplane  Y body  axis  component 
of  ai rspeed 

r 

0 

8 

5 

WBAS 

1/t 

Airplane  Z body  axis  component 
of  ai rspeed 

8 

6 

ALPHA 

rad 

Angle-of-attack 

1 1 

VJ 

8 

7 

BETA 

rad 

Sideslip  angle 

TT 

8 

8 

DUM1 

Dummy  variable,  this  location 

.. 

unused 

0 

8 

9 

DUM2 

Dummy  variable,  this  location 
unused 

8 

10 

V 5 

1/t 

Ship  speed 

D 

8 

11 

PSI 5 

rad 

Ship  heading 

n 

8 

12 

X1DWS 

1/t 

Earth  axis  component  of  ship's 

L 

relative  wind 

0 

8 

13 

Y1DWS 

1/t 

Earth  axis  component  of  ship's 
relative  wind 

8 

14 

VAMB 

1/t 

Ambient  wind  speed 

1 ^ 

8 

15 

PS IAMB 

rad 

Ambient  wind  heading 

1 

8 

16 

X1DWA 

1/t 

Earth  axis  component  of  ambient 
wind 

8 

17 

Y1DWA 

1/t 

Earth  axis  component  of  ambient 

i 

wind 

1 L 
■ 

8 

18 

Z1DWA 

1/t 

Earth  axis  component  of  ambient 

wind 

G 


~TF~ 


Cont’d 

R ARRAY 

SEGMENT  ELEMENT 

VARIABLE 

NAME 

UNITS 

i 

1 

1 

DESCRIPTION  1 

8 

19 

VW0D 

1/t 

Wind-over-deck  speed  X 

8 

20 

PSIW0D 

rad 

Wind-over-deck  heading  -- 

8 

21 

PSI3 

rad 

Orientation  of  ship  wing  axis; 
angle  of  rotation  about  the  Z 
earth  axis  to  point  ship  X wind 
axis  into  the  wind-over-deck 

8 

22 

PSISP 

rad 

= PSI3  - PSI5 

8 

23 

PSIS 

rad 

= -PSIP;  this  is  the  angle 

4»cS  in  the  airwake  math  model 
of  Reference  12 

l J 

8 

24 

DX45 

1 

Earth  axis  distance  between 
ai rplane  and  ship 

8 

25 

DY45 

1 

Earth  axis  distance  between 
airplane  and  ship 

»» 

8 

26 

DZ45 

1 

Earth  axis  distance  between 
airplane  and  ship 

8 

27 

XM 

* 

1 

Airplane  position  in  space  with 

8 

28 

YM 

1 

respect  to  the  ship  wind  axis  1, 

system  (i.e.,  the  airwake  model 

8 

29 

ZM 

1 . 

reference  point) 

8 

30 

PSIA 

rad  ' 

1 

Euler  type  angles  relating  the 

8 

31 

THETAA 

rad 

> airplane  to  the  ship  wind  axis 
system.  It  is  assumed  that  the 

8 

32 

PHIA 

rad  . 

ship  wind  axis  does  not  pitch 
and/or  roll . 

D 

8 

33 

SIGMAVX 

1/t  ' 

Standard  deviations  of  the 

8 

34 

SIGMAVY 

1/t 

> airwake  in  the  ship  wind  axis 
system 

8 

35 

SIGMAVZ 

vt. 

8 

36 

VXAW 

1/t' 

n 

Mean  components  of  the  ship 

8 

37 

VYAW 

1/t 

f 

> airwake  in  the  ship  wind  axis 
system 

8 

38 

VZAW 

1/t  J 

u 

J 


> 

-*• 



— - TT- - , 

!i 

Cont'd 

A 

R 

SEGMENT 

ARRAY 

ELEMENT 

VARIABLE 

NAME 

UNITS 

DESCRIPTION 

r 

0 

8 

39 

UWIND 

1/t 

Total  wind  component  along  air- 
plane X body  axis; 

UWIND  = UMEAN  + UAW 

0 

8 

40 

UMEAN 

1/t 

Ambient  wind  component  along 
airplane  X body  axis 

8 

41 

UR 

Random  wind  component;  unused 

Ml 

8 

42 

UAW 

1/t 

Mean  component  of  airwake 
disturbance  along  airplane  X 
body  axis 

0 

8 

43 

VWIND 

1/t 

Total  wind  component  along 
ai rplane  Y body  axis ; 

VWIND  = VMEAN  + VAW 

: 

11 

8 

44 

VMEAN 

1/t 

Ambient  wind  component  along 
airplane  Y body  axis 

8 

45 

VR 

Random  wind  component;  unused 

u 

0 

8 

46 

VAW 

1/t 

Mean  component  of  airwake 
disturbance  along  airplane  Y 
body  axis 

i 

0 

8 

47 

WWIND 

1/t 

Total  wind  component  along 
airplane  Z body  axis; 

WWIND  = WMEAN  + WAW 

1 B 

8 

48 

WMEAN 

1/t 

Ambient  wind  component  along 
airplane  Z body  axis 

8 

49 

MR 

Random  wind  component;  unused 

J 

8 

50 

WAW 

1/t 

Mean  component  of  airwake 
disturbance  along  airplane  Z 
body  axis 

u 

0 

8 

51 

WNSF 

rad/t 

Natural  frequency  of  second 
order  shaping  filter  used  in 
airwake  model 

I 

0 

8 

52 

PSISS 

sin  (PSIS) 

H 

E-17 

1 

— — - 

■ 

^ 

M^f — . 

_ 

— --r-  , 

1 

Cont ' d 


R ARRAY  VARIABLE 

SEGMENT  ELEMENT  NAME 


UNITS 


DESCRIPTION 


8 

53 

PSISC 

8 

54 

PSIST 

8 

55 

TS1 

8 

56 

TS2 

8 

57 

TS3 

8 

58 

TS4 

8 

59 

TS5 

8 

60 

TR1 

8 

61 

DXM145 

8 

62 

DXM245 

8 

63 

DYM145 

8 

64 

DYM245 

8 

65 

ZGSP 

8 

66 

ZBIASP 

cos  (PSIS) 
tan  (PSIS) 


Temporary  storage  for  preliminary 
calculations  involving  wind-over- 
^ deck  and  airwake  shaping  functions 


I Scale  factor  on  SIGMAVX, 

SIGMAVY,  and  SIGMAVZ 

= cos  (PSISP)*  cox  (PSI5)  - 
sin  (PSIP)*  sin  (PSI5) 

= cos  (PSIP)*  sin  (PS 1 5)  + 
sin  (PSIP)*  cos  ( PS 15) 

= -sin  (PSIP)*  cos  (PS 1 5)  - 
cox  (PSIP)*  sin  (PSI5) 

= cos  (PSIP)*  cos  (PSI5)  - 

sin  (PSIP)*  sin  (PSI5) 

1 Commanded  Z^p 

1 Bias  term 


E- 18 
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R ARRAY 

SEGMENT  ELEMENT 

VARIABLE 

NAME 

UNITS 

DESCRIPTION 

9 

1 

THETAC 

rad 

Coimanded  airplane  pitch 
attitude 

9 

2 

DFRPM0L 

rpm 

Open  loop  rpm  command 

9 

3 

DCDN0L 

rad 

Open  loop  nozzle  command 

9 

4 

DE0L 

units 

Open  loop  pitch  control 
command 

9 

5 

GAMMAGS 

deg 

Commanded  glide  slope  angle 

9 

6 

RANGE 

1 

- *$  - xA/p 

9 

7 

VAIRC 

1/t 

Commanded  airspeed 

9 

8 

DA0L 

units 

Open  loop  roll  control  command 

9 

9 

DR0L 

units 

Open  loop  yaw  control  command 

I 

J 


R 

SEGMENT 

ARRAY 

ELEMENT 

VARIABLE 

NAME 

UNITS 

DESCRIPTION  j 

1: 

1 

RC 

1 

Commended  rarge;  \ 

range  £ X^  - X^/p 

11 

2 

RCP 

1/t 

Commanded  range  tine  rate-of- 

change  ** 

11 

3 

XAPC0M 

I 

Airplane  commanded  position 
along  X earth  axis 

11 

4 

ERRX 

1 

= X^/p  - XAPC0M;  longitudinal 
position  error 

11 

5 

ERRU 

1/t 

A'rplane  body  axis  component  of 
error  rate 

LJ 

11 

6 

ERRZ 

1 

= - ZGSP;  vertical  position 

error 

11 

7 

EKRW 

1/t 

\ • 

Airplane  body  axis  component  of 

error  rate  r- 

11 

6 

ETHETA 

rad 

= THETA  - THETAC;  pitch  altitude 
err  o^ 

11 

9 

FRY 

1 

- YC0M  - v^/p;  lateral  position 
error 

11 

10 

ERPHI 

rad 

= PHIC0M  - PHI;  roll  altitude 
error 

11 

11 

ERPSI 

rad 

= PSIC0M  - PS I;  leading  error 

11 

12 

YC0M 

1 

Airplane  commanded  Y earth 
axis  position 

11 

13 

PHIC0M 

rad 

1 J 

Airplane  roll  cititude 

command  p 

11 

14 

PSIC0M 

rad 

Airplane  heading  command 

0 

D 

II 


"CL. 


I 

I 

l 

r; 

l! 


11 


L 


R 
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VARIABLE 
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UNITS 
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12 

1 

DE 

units 

Pitch  control  displacer, ert 

12 

2 

DT 

units 

Thrust,  control  displacement 

12 

3 

DCDN 

rad 

Nozzle  position 

12 

4 

DED 

1/t 

Pitch  control  rate 

12 

5 

DTD 

1/t 

Thrust  control  rate 

12 

6 

DCDND 

rad/t 

Nozzle  angle  rate 

12 

7 

DA 

units 

Roll  control  displacement 
limited  output  value 

12 

8 

DR 

units 

Yaw  control  displacement 
limited  output  value 

12 

9 

DAP 

units 

Roll  control  displacement, 
input  to  limiter 

12 

10 

DRP 

units 

Yaw  control  displacement, 
input  to  limiter 

12 

11 

DAI 

units 

Intermediate  variable  used  in 
calculation  of  roll  control 

12 

12 

DA  2 

units 

Intermediate  variable  used  in 
calculation  of  roll  control 

12 

13 

DR  1 

units 

Intermediate  variable  used  in 
calculation  of  yaw  control 

12 

14 

DR2 

units 

Inter-mediate  variable  used  in 
calculation  of  yaw  control 

12 

15 

DA  ID 

1/t 

Intermediate  rate  variable  used 
in  calculation  of  roll  control 

12 

16 

DA2D 

1/t 

Intermediate  rate  variable  used 
in  calculation  cf  roll  control 

12 

17 

DR1D 

1/t 

Intermediate  rate  variable  used 
in  calculation  of  yaw  control 

12 

18 

DR2D 

1/t 

Intermediate  rate  variable  used 
in  calculation  of  yaw  control 

* 


E- 23 
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13 

1-6 

TABV1 

1/t 

Independent  variable  table 
airspeeds 

of 

14 

1-10 

TABV2 

t 

Independent  variable  table 
time 

of 

E-24 


R ARRAY 


VARIABLE 


UNITS 


DFSCRIPTION 


SIGDAP 


SIGDRP 


SI.GDA 


SI3DR 


GAIN DA 


GAINDR 


units 


units 


uni  ts 


units 


Standard  deviation  of  roll 
control  displacement  at  input 
to  limiter 

Standard  deviation  of  y3w 
control  displacement  at  input 
to  limiter 

Standard  deviation  of  roll 
control  displacement  after 
1 imi ter 

Standard  deviation  of  yaw 
control  displacement  after 
limiter 

Gain  on  random  component  of 
roll  control,  from  describing 
function  representation  of 
limiter 

Gain  on  random  component  of 
yaw  control,  from  describing 
function  representation  of 
1 imi ter 


R ARRAY 
SEGMENT  ELEMENT 


VARIABLE 

NAME 


UNITS 


DESCRIPTION 


KXS 

1 

Gains,  damping  ratios,  and 

1 

natural  frequency  for  the 

1 

ZETAXS 

transfer  functions  which 

describe  ship  surge,  heave, 

7 

WNXS 

ra  d/t 

pitch,  sway,  roll,  and  yaw 

1 

perturbations 

i 

ZETAYS 

| 

WNYS 

rad/t 

i 

KRS 

rad 

1 

ZETAPS 

1 

WNRS 

rad/t 

dl 

KFS 

rad 

3 

ZETAPS 

. 

WNPS 

rad/t 

11 

TAUPRT 

t 

Time  constant  of  first  o^der 
lac.  Output  of  this  filter 
modifies  pilots  altitude 
command. 

0 

D 

0 

0 

E-26 

H 

f 


I 
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UNITS 

DESCRIPTION 

11 

18 

None 

N/A 

N/A 

Open  segment 

Q 

19 

None 

N/A 

N/A 

Open  segment 

0 

20 

1-10 

TABTC 

rad 

Dependent  variable  table  of 
pitch  altitude  commands 

L 

21 

1-10 

TRPM0L 

rpm 

Dependent  variable  table  of 
rpm  open  loop  inputs 

1 

L_ 

22 

1-10 

TCDN0L 

deg 

Dependent  variable  table  of 
nozzle  position  open  loop 
inputs 

D 

23 

1-10 

TDEPL 

units 

Dependent  variable  table  of 
pitch  control  open  loop  inputs 

f] 

il 

24 

1-10 

TABR 

Dependent  variable  table  of 
range  commands;  net  currently 
in  use 

c 

25 

1-10 

TABRD0T 

Dependent  variable  table  of 
rate-of-closure;  not  currently 
in  use 

[] 

0 

27 

1-6 

TX'J 

Dependent  variable  tables  of 
longitudinal  dimensional 
derivatives.  Units  are  consistent 
with  output  values 

28 

1-6 

TXW 

D 

29 

1-6 

TXQ 

n 

30 

1-6 

TXR°M 

u 

31 

1-6 

TXCDN 

n 

32 

1-6 

TXPIGV 

U 

33 

1-6 

TXDE 

C 

34 

1-6 

TZU 

35 

1-6 

TZW 

I 

1 1 

A 

R 

ARRAY 

VARIABLE 

u 

SEGMENT 

ELEMENT 

NAME 

UNITS 

DESCRIPTION 

rj 

50 

varies 

F 

Overall  system  matrix 

1 * 

51 

varies 

P 

Covariance  matrix 

[ 

52 

varies 

PD0T 

Covariance  matrix  time  rate-of- 
change;  PD0T=F*P  + P*FT  + GM*QM*GMT 

53 

varies 

GM 

Matrix  relating  external  white 
noise  disturbances  tc  state 
equations 

5^ 

varies 

QM 

Covariance  matrix  of  external 
white  noise  sources 

: L 

55 

1-50 

SIGMA 

varies 

Array  of  standard  deviations 

56 

1-50 

SIGMAXY 

varies 

Array  of  cross  correlations  with 

a given  element  (used  when 
plotting  ellipses);  see  subroutine 


miscal 


Li 

57 

1-75 

vxx 

1/t 

Table  of  mean  components  of 
airwake  along  x ship  wind  axis 

D 

58 

1-75 

VYY 

1/t 

Table  of  mean  components  of 
a invoke  along  y ship  wind  axis 

0 

59 

1-75 

VZZ 

1/t 

Table  of  mean  component  of 
airwake  a’ong  z ship  wind  axis 

60 

1-75 

S VXX 

1/t 

Table  of  standard  deviation 
components  of  airwake  along 

X ship  wind  axis 

0 

61 

1-75 

SVYY 

1/t 

Table  of  standard  deviation 
components  of  airwake  3long  Y 
ship  wind  axis 

n 

62 

1-75 

SVZZ 

1/t 

Table  of  standard  deviation 

components  along  Z ship  wind 
axis 


I 

I 

0 

1 


I 


R 

SEGMENT 

ARRA'i 

ELEMENT 

VARIABLE 

NAME 

UNITS 

description 

63 

1-6 

TKDEU 

OeDendent  variable  tables  of 
pilot  feedback  ga-ns.  Units 
are  consistent  with  oreviousl 
defined  output  values. 

64 

1-6 

TKDEW 

65 

1-6 

TKDEQ 

66 

1-6 

TKDET 

I 

67 

1-6 

TKDEX 

l 

68 

1-6 

TKDEZ 

69 

1-6 

TKDEN 

70 

1-6 

TKDEDE 

71 

1-6 

TKDEDT 

72 

1-6 

TKDTU 

73 

1-6 

TKDTW 

74 

1-6 

TKDTQ 

75 

1-6 

TKOTT 

76 

1-6 

TKDTX 

77 

1-6 

TKDTZ 

78 

1-6 

TKDTN 

79 

1-6 

TKDTDE 

80 

1-6 

TKDTDT 

81 

varies 

Fl 

Aircraft  system  matrix 

82 

varies 

G1 

Aircraft  control  matrix 

83 

varies 

GAMMA 

Matrix  relating  colored  noise 
inputs  to  aircraft 

T 


11 


1 

1 

I 

1 


\ 


fl 

0 

0 

fl 

1 
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R ARRAY 

SEGMENT  ELEMENT 

VARIABLE 

NAME 

UNITS 

DESCRIPTION 

84 

varies 

A 

Colored  noise  shaping  matrix 

85 

varies 

B 

Matrix  relating  white  noise 
inputs  to  shaping  filters 

86 

varies 

Ml 

Not  currently  used 

87 

varies 

M2 

Not  currently  used 

88 

varies 

H 

Not  currently  used 

89 

varies 

D 

Measurement  matrix 

90 

varies 

C 

Feedback  matrix  operating  on 
states  in  control  rate  equation 

91 

varies 

RL 

Feedback  matrix  operating  on 
controls  in  control  rate 
equation 

92 

varies 

W 

Covariance  matrix  of  white  noise 
external  disturbances 

93 

varies 

VU 

Covariance  matrix  of  white 
neuromuscular  noise 

94 

varies 

VY 

Covariance  matrix  of  white 
measurement  noise 

95 

varies 

VYI 

Inverse  of  matrix  VY 

96 

varies 

S 

Ship  motion  system  matrix 

97 

varies 

FE 

System  matrix  for  estimator 
covariance  calculation 

98 

va  ri  es 

QE 

Covariance  matrix  for  driving 
estimator  covariance  calculation 

99 

varies 

PE 

Covariance  matrix  for  estimator 
error 

100 

varies 

PED 

Time  rate-of-change  of  estimator 
error  covariance  matrix 

101 

varies 

KG 

Kalman  gain  matrix 
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R 

SEGMENT 

ARRAY 

ELEMENT 

VARIABLE 

NAME 

UNITS 

DESCRIPTION 

102 

va  ri  es 

FXHAT 

System  matrix  for  estimated 
state  covariance  calculations 

103 

varies 

PXHAT 

Estimated  states  covariance 
matrix 

104 

varies 

PXHATD 

Time  rate-of -change  of  estimated 
states  covariance  matrix 

105 

varies 

PX 

State  covariance  matrix 

106 

varies 

CS 

A sub  matrix  in  F related  to  the 
modified  pilot  altitude  command 
due  to  ship  motion 

R ARRAY 

SEGMENT  ELEMENT 


VARIABLE 

NAME 


UNITS 


DESCRIPTION 


107 

1 

DEUR 

107 

2 

DEWR 

107 

3 

DEVX 

107 

4 

DEVXD 

107 

5 

DEVZ 

107 

6 

DEVZD 

107 

7 

DTUR 

107 

8 

DTWR 

107 

9 

DTVX 

107 

10 

DTVXD 

107 

11 

DTVZ 

107 

12 

DTVZD 

Y These  additional  feedback  gains 
in  the  optimal  pilot  model 
operate  on  the  estimated  colored 
noise  disturbances;  results  of 
table  look-up 

1 

T 

1 

T 

t 

1 

I 

1 

t 

1 

i 

1 

1 

1 

1 

1 

t 

1 


t 

1 


I 

IS 


r 


R ARRAY 

SEGMENT  ELEMENT 


VARIABLE 

NAME 


UNITS 


DESCRIPTION 


108 

1-6 

TDEUR 

109 

1-6 

TDEWR 

no 

1-6 

TDEVX 

111 

1-6 

TDEVXD 

112 

1-6 

TDEVZ 

113 

1-6 

TDEVZD 

114 

1-6 

TDTUR 

115 

1-6 

TDTWR 

116 

1-6 

TDTVX 

117 

1-6 

TDTVXD 

118 

1-6 

TDTVZ 

119 

1-6 

TDTVZD 

120 

1-6 

TYV 

121 

1-6 

TYP 

122 

1-6 

TYR 

123 

1-6 

TYDA 

124 

1-6 

TYDR 

125 

1-6 

TLV 

126 

1-6 

TLP 

127 

1-6 

TLR 

Dependent  variable  tables  for 
pilot  feedback  gains.  Units 
are  consistent  with  the  output 
val ues . 


Dependent  variable  tables  for 
lateral /directional  dimensional 
derivatives.  Units  are  consistent 
with  output  values. 


] 


3 


1 

] 

3 


is 


n 

c 

E 

» t 
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3 
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ARRAY 

ELEMENT 

VARIABLE 

NAME 

UNITS 

128 

1-6 

TLDA 

129 

1-6 

TLDR 

130 

1-6 

TNV 

131 

1-6 

TNP 

132 

1-6 

TNR 

133 

1-6 

TNDA 

134 

1-6 

TNDR 

135 

1-10 

TABUB 

1/t 

136 

1-10 

TABWB 

1/t 

137 

1-10 

TXAPP 

1 

138 

1-10 

TZAPP 

1 

139 

1-10 

TXSP 

1 

140 

1-4 

TDAIN 

units 

141 

1-4 

mA0UT 

units 

142 

1-4 

TDR1N 

uni  ts 

143 

1-4 

TDR0UT 

units 

Dependent  variable  tables  of 
longitudinal  parameters  required 
in  a lateral /directional  run. 


Independent  variable  table  for 
defining  roll  control  limiter. 

Dependent  variable  table  for 
defining  roll  control  limiter. 

Independent  variable  table  for 
defining  yaw  control  limiter. 

Dependent  variable  table  for 
defining  yaw  control  limiter. 
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F.O  LIST  OF  SYMBOLS 

Throughout  the  text  several  equations  appear  in  vector-matrix  format. 

In  these  equations  lower  case  letters  represent  vectors  and  upper  case 
letters  are  matrices.  For  scalar  quantities,  upper  case  letters  are  nominal 
or  mean  values  and  lower  case  letters  represent  perturbation  values. 

The  test  clearly  distinguishes  between  variables,  as  to  whether  they  are 
matrices,  vectors,  or  scalars,  when  they  are  used. 


E { } 


Dimensional  stability  derivative 
(includes  mass  or  inertia); 

C = X,  Y,  Z,  L,  M,  and  N; 
i = u,  v,  w,  p,  q,  r,  and  <5j . 

Expected  value  operator. 

Gravitational  constant. 


y ' y 1 7 1 

A/P’  A/P’  Lk/ P 


Xj,  Yj,  Zj 


W.O.D. 


Laplace  transform  variable. 

Earth  axis  coordinates  of  the 
airplane. 

Earth  axis  coordinates  of  the 
ship. 

Wind  over  deck 
Thrust  to  weight  ratio 
Deflection  of  j—  controller. 
AV-8  pitch  control  position. 
AV-8  throttle  position. 

AV-8  nozzle  angle. 

AV-8  roll  control  position 
AV-8  yaw  control  position 
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o 


Standard  deviation. 


Euler  angles. 

Time  derivative  of  ( ) . 
Mean  value  of  ( ). 

Expected  value  of  ( ) . 

Denotes  reference  value  or 
condition. 

Gust. 
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